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Abstract. We discuss the constraints set on galaxy evolution by a variety of data from deep extragalactic surveys 
performed in the mid-IR and far-IR with the Infrared Space Observatory and with millimetric telescopes at 
longer wavelengths. These observations indicate extremely high rates of evolution for IR galaxies, exceeding those 
measured for galaxies at other wavelengths and comparable or larger than the rates observed for quasars. We also 
match the modelled integrated emission by IR galaxies at any redshifts with the observed spectral intensity of 
the extragalactic IR background (CIRB), as a further constraint. The multi-wavelength statistics on IR galaxies 
can be reconciled with each other by assuming for the bulk of the population spectral energy distributions (SED) 
as typical for starbursts, which we take as an indication that stellar (rather than AGN, see also Fadda et al. 
2001) activity powers IR emission by faint galaxies. According to our model and following the analysis of Elbaz 
et al. (2001), the deep ISO surveys at 15 nm may have already resolved more than 50% of the bolometric CIRB 
intensity: the faint ISO 15 /im source samples, relatively easy to identify in deep optical images (Aussel et al. 
1999), can then allow to investigate the origin of the CIRB background. From our fits to the observed optical-IR 
SEDs, these objects appear to mostly involve massive galaxies hosting luminous starbursts {SFR ~ 100 Mq/jt). 
The evolutionary scheme we infer from these data considers a bimodal star formation (SF), including a phase of 
long-lived quiescent SF, and enhanced SF taking place during transient events recurrently triggered by interactions 
and merging. We interpret the strong observed evolution as an increase with z of the rate of interactions between 
galaxies {density evolution) and an increase of their IR luminosity due to the more abundant fuel available in the 
past (luminosity evolution): both factors enhance the probability to detect a galaxy during the "active" phase at 
higher z. Very schematically, we associate the origin of the bulk of the optical/NIR background to the quiescent 
evolution, while the CIRB is interpreted as mostly due the dusty starburst phase. The latter possibly leads to the 
formation of galaxy spheroids, when the dynamical events triggering the starburst re-distribute already present 
stellar populations. The large energy contents in the CIRB and optical backgrounds are not easily explained, 
considering the moderate efficiency of energy generation by stars: a top-heavy stellar IMF associated with the 
starburst phase (and compared with a more standard IMF during the quiescent SF) would alleviate the problem. 
The evolution of the IR emissivity of galaxies from the present time to z ~ 1 is so strong that the combined set 
of constraints by the observed z-distributions and the CIRB spectrum impose it to turn-over at z > 1: scenarios 
in which a dominant fraction of stellar formation occurs at very high-z are not supported by our analysis. 
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1. Introduction 

High-redshift galaxies and the generation of stars dur- 
ing the past cosmic history are most usually investi- 
gated by observing the stellar integrated emissions in the 
UV/optical/near-IR. These studies, based on a variety of 
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selection techniques and exploiting very large telescopes 
on ground and efRcient photon detectors, achieved ex- 
traordinary successes during the last ten years in dis- 
covering large numbers of high-z galaxy candidates, most 
of which were later confirmed by high-sensitivity optical 
spectroscopy (Madau et al. 1996; Steidel et al. 1999). The 
outcomes of these unbiased surveys are galaxies charac- 



2 



Franccschini A. et al.: Galaxy Evolution at Long Wavelengths 



terized by typically moderate luminosities and a modest 
activity of star- formation (few M0/yr on average). 

However, that this view could be to some extent incom- 
plete is illustrated by the fact that biased optical surveys 
emphasize a quite more active universe, as revealed by the 
existence of heavily metal-enriched environments around 
quasars and active galaxies at any redshifts (Omont et al. 
1996; Padovani & Matteucci 1993; Pranceschini & Gratton 
1997); by the presence of populations of massive elliptical 
galaxies up to 2: > 1, with dynamically relaxed profiles and 
complete exhaustion of the ISM (e.g. Moriondo, Cimatti, 
Daddi 2000; Rodighicro, Franccschini, Fasano 2001), ex- 
pected to originate from violent star bursts; and by the 
intense activity of massive stars required to explain the 
metal-polluted hot plasmas present in galaxy clusters and 
groups (Mushotzky, Loewenstein, 1997). What is essen- 
tially missing from optical observations is the evidence of 
cosmic sites where active transformations of baryons are 
taking place at rates high enough to explain the above 
findings, among others. 

Hints on such a possible missing link between the ac- 
tive and quiescent universe come from inspection of the lo- 
cal universe. The IRAS long- wavelength surveys, in partic- 
ular, have revealed that in a small fraction of local massive 
galaxies (the so-called luminous [LIRG] and very luminous 
[ULIRG] infrared galaxies) star-formation is taking place 
at very high rates (> 100 Mg/yr) (Sanders et al. 1988; 
Kormendy & Sanders 1992). Interesting to note, the red- 
dened optical spectra of these objects do not contain mani- 
fest signatures of the dramatic phenomena revealed by the 
far-IR observations (Poggianti and Wu 1998). All this em- 
phasizes the role of extinction by dust, which is present 
wherever stars are formed, but is increasingly important in 
the most luminous objects. This also illustrates the power 
of extending the selection waveband for cosmological sur- 
veys from the optical (tracing the stellar-dominated emis- 
sion) to the IR where ISM-dominated emission is observ- 
able. 

It is only during the last five years that new powerful 
instrumentation has allowed to start a systematic explo- 
ration of the distant universe at long wavelengths. Three 
major developments have allowed this. Firstly, the dis- 
covery in the COBE all-sky maps of a bright isotropic 
background in the far-IR/sub-mm, of likely extragalactic 
origin (CIRB) and interpreted as the integrated emission 
by dust present in distant and primeval galaxies (Puget 
et al. 1996; Hauser et al. 1998). 

The second important fact was the start of operation 
of the bolometer array SCUBA on the 15m sub-millimetric 
telescope JCMT, able to resolve a substantial fraction of 
the CIRB background at long wavelengths into a popu- 
lation of very luminous IR galaxies at z ~ 1 or larger 
(Small et al. 1997; Hughes et al. 1998; Barger et al. 1998; 
Blain et al. 1999). A new powerful bolomctric imaging 
camera (MAMBO) has also recently become operative on 
the IRAM 30m telescope, and started to provide clean 
deep images of the distant universe at 1300 fxm (Bertoldi 
et al. 2000, 2001). 



Finally, the Infrared Space Observatory (ISO) allowed 
for the first time to perform sensitive surveys of distant IR 
sources in the mid- and far-IR (Elbaz et al. 1999; Puget et 
al. 1999) and to characterize in detail the evolution of the 
IR emissivity of galaxies up to redshift 2; ~ 1 and above. 

The present paper is devoted to the analysis of a large 
dataset including number counts, redshift distributions, 
and luminosity functions for faint IR sources selected be- 
tween Aeff — 10 and ~ 1000 //m. Wc emphasize in our 
analysis deep survey data from the ISO mission, partic- 
ularly in the mid-IR where the number count statistics 
arc robust (Elbaz et al. 1999; Altieri, Metcalfe and Kneib 
1999; Aussel et al. 1999; Oliver et al. 1997). 

Our approach is different from those of previous mod- 
els fitting the IR galaxy counts. In most cases (e.g. 
Devriendt et al. 1999, Roche and Eales 1999; Pearson & 
Rowan- Robinson 1996; Rowan- Robinson 2001; Xu et al. 
2001) attempts have been made to provide combined de- 
scriptions of the IR and optical-UV data on faint galax- 
ies. This approach could produce even misleading results 
whenever the optical data would constrain the global so- 
lution to give very poor fits of the IR data. Only an ex- 
tremely detailed description of the complex relationship 
between optical and IR emissions could provide meaning- 
ful results at some stages. In our view, published models 
(e.g. Guiderdoni et al. 1998, including sophisticated mod- 
elling of the formation of structures) illustrate more the 
inconsistencies emerging when comparing optical and IR 
statistics on faint sources than the benefits of a combined 
analysis. 

Another way to see the problem is to consider in 
some details the optical spectral properties of luminous 
and ultra- luminous IR galaxies. Poggianti & Wu (2000), 
Poggianti, Bressan & Franceschini (2001) and Rigopoulou 
et al. (2000) have studied rest-frame optical spectra for 
both local and high-redshift objects, and consistently 
found that ~ 70% — 80% of the energy emitted by young 
stars and reprocessed in the far-IR leaves no traces in the 
optical spectrum (even after correction for dust extinc- 
tion) , hence can only be accounted for by long- wavelength 
observations. Altogether, the ratio of IR to optical emis- 
sions is very broadly distributed and no clear empirical, 
nor physical, relationships have yet been established be- 
tween the two. 

Consequently, we have chosen to confine our analysis 
to data at long-wavelengths (10 to 1000 /um) and to search 
for simple parametrizations of the evolution of galaxies in 
the IR, as an attempt to provide guidelines for future phys- 
ical models of galaxy activity and its evolution. Only at 
the end we will relate these results with optical data on 
faint galaxies, by comparing the observed integrated emis- 
sions in the forms of the optical and CIRB backgrounds. 
In spite of the schematicity of such an approach, our re- 
sults already contain new critical information on galaxy 
formation and evolution. 

The paper is organized as follows. We set first the 
general framework by discussing in Sect. 2 the informa- 
tion contained in the CIRB spectral intensity. Indeed, the 
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availability of the CIRB measurements, providing a solid 
constraint on the integrated IR emission of galaxies at any 
epochs, is a rather unique feature of the IR domain, mak- 
ing it of extreme interest for studies of galaxy evolution. 
We mention here some recent measurements of the cosmic 
opacity at very high (TeV) energies, allowing to set rele- 
vant constraints at wavelengths where the CIRB is not di- 
rectly measurable. Sect. 3 is devoted to summarize results 
from the most relevant survey projects with ISO, while in 
Sect. 4 some relevant data obtained with millimetric tele- 
scopes are summarized. Sect. 5 illustrates our attempt to 
reproduce the multi- wavelength data with simple prescrip- 
tions. Our present understanding of the physical nature of 
the IR source populations is discussed in Sect. 6, together 
with a simple physical interpretation of their previously 
described evolution. Sect. 7 contains a discussion of the 
global properties of high-redshift IR galaxies, like the evo- 
lutionary SFR density, and of the energy constraints set 
by the observed IR and optical backgrounds. Our conclu- 
sions are summarized in Sect. 8. 

For consistency with previous analyses, we adopt for 
Ho the value of 50 Km/sec/Mpc (note that this choice 
has no impact on our inferred evolution properties, since 
the dependences on Hq of the luminosity functions and 
number counts cancel out). In the following we indicate 
with the symbol L12 the luminosity i'L{i') calculated at 
A = 12/im and expressed in solar units. The same termi- 
nology is used for L at other wavelengths. The symbol S12 
indicates the monochromatic flux (in Jy) at 12 fXTn (and 
similarly for other wavelengths). 

2. SETTING THE FRAMEWORK: THE COSMIC 
INFRARED BACKGROUND 

Cosmic background radiations provide a fundamental 
channel of information on high-redshift sources, partic- 
ularly when, for technological limitations, observations at 
faint flux levels in a given waveband are not possible (as it 
is largely the case in the IR/sub-mm domain). We briefly 
review in this Section the observational status about the 
recently discovered cosmological background at IR and 
sub-millimetric wavelengths (CIRB), providing important 
constraints on galaxy formation and evolution. 

The discovery of the CIRB - anticipated by a model- 
listic prediction by Franceschini et al. (1994), and made 
possible by the NASA's COBE mission - was viewed as 
the Erst chance to determine, or at least constrain, the 
integrated emission of distant galaxies (Puget et al. 1996; 
Guiderdoni et al. 1997; Hauser et al. 1998; Fixsen et al. 
1998). For comparison, extragalactic backgrounds at other 
wavelengths appear to contain only moderate contribu- 
tions by distant galaxies: the Radio, X-ray and 7-ray back- 
grounds, apparently dominated by distant quasars and 
AGNs (Giacconi et al. 2001; Tozzi et al. 2001), and the 
Cosmic Microwave Background including photons gener- 
ated at z ~ 1500. Also, direct measurements of the optical- 
UV backgrounds are hampered by the intense starlight 



reflected by high latitude "cirrus" dust and Zodiacal- 
reflected Sun-light. 

2.1. Observational status about the CIRB 

In spite of the presence of bright foregrounds (Zodiacal 
and Interplanetary dust emission. Galactic Starlight, high- 
latitude "cirrus" emission), there are two relatively clean 
spectral windows in the IR where these summed emissions 
produce two minima: the near-IR (3-4 /um) and the sub- 
mm (100-500 /Ltm) cosmological windows. Redshifted pho- 
tons by the two most prominent galaxy emission features, 
the stellar photospheric peak at A 1 /im and the one at 
A 100 /im due to dust re-radiation, are here observable 
in principle. 

Particularly favourable for the detection of an extra- 
galactic signal turned out to be the longer-wavelength 
channel. By exploiting the different spatial dependen- 
cies of the various dust components and the observed 
correlations with appropriate dust tracers like the neu- 
tral and ionized hydrogen (throug the HI 21 cm and H^ 
lines), Puget et al. (1996) have identified in the all-sky 
FIRAS/COBE maps an isotropic signal with an intensity 
following the law vB^ ~ 3.4 x 10^^(400 ^J,m/Xf W m^^ 
sr^^ in the 400-1000 fim interval. 

This tentative detection has been confirmed by vari- 
ous other groups with independent analyses of data from 
FIRAS on COBE (e.g. Fixsen et al. 1998), as well as from 
the DIRBE experiment in two broad-band channels at 
A = 140 and 240^r7i (Hauser et al. 1998). The most re- 
cent results from DIRBE (Hauser et al. 1998; Lagache et 
al. 1999; Finkbeiner, et al. 2000) and FIRAS (Fixsen et 
al. 1998) are reported in Fig. ^ 

Finkbeiner, Davies & Schlegel (2000), after a very del- 
icate subtraction of the far dominant Galactic and IPD 
foregrounds, found an isotropic signal at 60 and 100 /xm 
with intensities at the level of ~ 30 10^^ W m^^ sr^^. 
This controversial result (see Puget & Lagache 2001 for a 
critical assessment) appears to conflict, in any case, with 
independent estimates based on observations of the cosmic 
high-energy opacity (see below). 

Recent analyses by Dwek & Arendt (1998) and 
Gorjian, Wright & Chary (2000) have claimed tentative 
detections in the near-IR cosmological window at 3.5 fj,m 
and in the J, H and K DIRBE bands, however with large 
uncertainties because of the very problematic evaluation 
of the Zodiacal scattered light. Because of this, CIRB es- 
timates particularly in J, H and K are to be taken more 
reliably as upper limits. 

2.2. Constraints from observations of the cosmic 
high-energy opacity 

No significant isotropic signals were detected at 4/im < 
A < 60^m, any cosmological flux being far dominated 
here by the Interplanetary dust (IPD) emission (to reduce 
it, missions to the outer Solar System would be needed). 
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Fig. 1. The Cosmic Infrared Background (CIRB) spectrum as measured by independent groups in the all-sky COBE 
maps (e.g. Hauser et al. 1998), compared with estimates of the optical extragalactic background based on ultradeep 
optical integrations by the HST in the HDF (Madau & Pozzetti 2000). The three lower datapoints in the far-IR are 
from a re-analysis of the DIRBE data by Lagache et al. (1999), the shaded areas from Fixsen et al. (1998) and Lagache 
et al. The two mid-IR points are the resolved fraction of the CIRB by the deep ISO surveys IGTES (Elbaz et al. 2001), 
while the dashed histograms are limits set by TeV cosmic opacity measurements (Sect. 2.2). The lower dashed line is 
the expected intensity based on the assumption that the IR emissivity of galaxies does not change with cosmic time. 
The thick line is the predicted CIRB spectrum of the presently discussed reference model for IR galaxy evolution. The 
dotted line marked CBR corresponds to the Cosmic Microwave Background spectrum. 



In this wavelength interval the CIRB energy density can 
be presently constrained with high-energy observations of 
Blazars, by measuring the optical depth at TeV energies 
due to the 7^7 interaction with the background CIRB 
photons (Stecker, de Jager & Salomon 1992). 

The absorption cross-section of 7-rays of energy E^ 
has a maximum for IR photons with energies obeing the 
condition: Cmax = 2{meC^Y /E^, or equivalently: Apoak — 
1.24 ± 0.6(£;^[TeV]) nm. The optical depth for a high- 
energy photon Eq travelling through a cosmic medium 
filled of low-energy photons with density p{z) from Zc to 
the present time is 

t{Eo,z,)= (1) 

c I dz— I dx— [ dii{\ + zf' ^'^j^ ^ a^^{E(i,i'[\ + z]) 
Jo dz Jq 2 Jq hv 



where (7 -y'-y IS the cross-section for photon-photon interac- 
tion. Coppi & Aharonian (1999) report the following ana- 
lytical approximation, good to better than 40%, to eq.(|l]): 

Applications of this concept have been possible when 
data from the Gamma Ray Observatory and X-ray space 
telescopes have been combined with observations at TeV 
energies by the Whipple, HEGRA and other Cherenkov 
observatories on Earth. Stanev & Franceschini (1998) have 
discussed upper limits on the CIRB with minimal a-priori 
guess on the CIRB spectrum, using HEGRA data for the 
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Blazar MKN 501 (z=0.034) during an outburst in 1997, 
on the assumption that the high-energy source spectrum 
is the flattest aUowed by the data. These hmits (dotted 
histogram in Fig. |l|) get quite close to the background 
flux already resolved by the ISO mid-IR deep surveys (see 
Sect. 5.4 below). 

More recently, Krawczynski et al. (1999) have com- 
bined the observations by Aharonian et al. (1999) of the 
MKN501 1997 outburst with X-ray data from RossiXTE 
and BeppoSAX, providing a simultaneous high-quality de- 
scription of the whole high-energy spectrum. These data 
are well fit by a Synchrotron Self Compton (SSC) model 
in which the TeV spectrum {i/ ^ 10^^ Hz) is produced by 
Inverse Compton of the hard X-ray spectrum {v ~ 10^* 
Hz): the combination of the two constrains the shape of 
the "primary" (i.e. before cosmic attenuation) spectrum 
at TeV energies. This is used to derive r^-y as a function 
of energy and, after eqs. (|l|) and (||), a constraint on the 
spectral intensity of the CIRB. The result is compatible 
with the limits by Stanev & Franceschini (1998, see also 
Renault et al. 2001) and allows a tentative, model depen- 
dent, estimate of the CIRB intensity in the interval from 
A = 20 to 40 /Ltm (see Fig.§). 

The observations of purely power-law Blazar spectra 
around ~ 1 TeV translate into a fairly robust upper 
limit of about 10~^ Watt/m^ sr at A ~ l/im shown in 
Fig. |l|. Substantially exceeding it, as sometimes suggested 
(Bernstein et al. 1998, Gorjian et al. 2000), would imply 
either very "ad hoc" 7— ray source spectra or new physics 
(Harwit, Proteroe & Bierman 1999). 

2.3. Constraints by CIRB and optical backgrounds on 
galaxy evolution 

Altogether, after years of active debate among various 
teams working on the COBE data, first about the ex- 
istence and later on the intensity and spectral shape of 
CIRB, there is now ample consensus, at least from 140 
to 500 /im where the CIRB spectrum is most reliably 
measured and where two completely independent datasets 
(FIRAS and DIRBE, with independent absolute calibra- 
tions) are available. The CIRB flux has in particular sta- 
bilized at values i^Ii, ~ 20 ± 5 and 1^ ~ 15 ± 5 10^^ 
Watt/m^/sr at A = 140 and 240 ^m. Modest differences 
in the calibration of FIRAS and DIRBE around 100 /zm 
have been reported (Hauser et al. 1998), but these do not 
affect the overall result. 

The measurement of the CIRB provides the global en- 
ergy density radiated in the IR by cosmic sources at any 
redshifts. Two concomitant facts - the very strong K- 
correction for galaxies in the far-IR/sub-mm due to the 
very steep and featureless dust spectra, and their robust- 
ness due to the modest dependence of dust equilibrium 
temperature T on the radiation field intensity - have sug- 
gested to use the CIRB spectrum to infer the evolution 
of the long-wavelength galaxy emissivity as a function of 
redshift (Gisper, Lagache & Puget 2000). Indeed, while 



the peak intensity at A = 100 to 200 /im constrains it 
at z < 1, the low foreground contamination at A > 200 
/im allows to set important constraints on the universal 
emissivity at z > 1. 

Between 100 and 1000 /im the observed integrated 
CIRB intensity turns out to be - (30 ± 5) 10"^ 
Watt/m^/sr. In addition to this measured part of the 
CIRB, one has to consider the presently un-measurable 
fraction resident between 100 and 10 /tni. Adopting mod- 
ellistic extrapolations as in Fig. |l|, consistent with the con- 
straints set by the cosmic opacity observations, the total 
energy density between 7 and 1000 rises to 

z//(z^)|fir 40 10-^ Watt/mVsr. (3) 

This flux is to be compared with the integrated bolo- 
metric emission by distant galaxies between 0.1 and 7 /tm 
(the "optical background"), for which we adopt the value 
given by Madau & Pozzetti (2000): 

z//(t/)|opt =i (17 ± 3) 10"^ Watt/mVsr. (4) 

This latter has been obtained from HST source counts be- 
tween 0.3 and 3 /tm down to the faintest detectable limits, 
by exploiting the number count convergence at magni- 
tudes mAB > 22. A significant upwards revision of this 
optical background suggested by Bernstein et al. [1998] to 
account for low surface brigtness emission by galaxies is 
not confirmed and would tend to conflict with measure- 
ments of the cosmic high-energy opacity. 

Already the directly measured part of the CIRB sets 
a relevant constraint on the evolution of cosmic sources, 
if we consider that for local galaxies only 30% on average 
of the bolometric flux is absorbed by dust and re-emitted 
in the far-IR. The CIRB's intensity exceeding tlie opti- 
cal background tells that galaxies in the past should have 
been much more "active" in the far-IR than in the optical, 
and very luminous in an absolute sense. A substantial frac- 
tion of the whole energy emitted by high-redshift galaxies 
should have been reprocessed by dust at long wavelengths. 

3. RESOLVING THE CIRB INTO SOURCES: 
DEEP SKY SURVEYS WITH THE 
INFRARED SPACE OBSERVATORY (ISO) 

The ISO Observatory (a 60cm cryogenic telescope oper- 
ated by ESA between 1995 and 1998) included two in- 
struments of cosmological interest: a mid-IR 32x32 ar- 
ray camera (ISOCAM), and a far-IR imaging photometer 
(ISOPHOT) with small 3x3 and 2x2 detector arrays from 
60 to 200 /tm. The main extragalactic results from the 30- 
month ISO mission have been summarized by Genzel & 
Cesarsky (2000). 

The improvement in sensitivity offered by ISO with 
respect to the previous IRAS surveys motivated to spend 
a relevant fraction of the observing time to perform a set 
of deep surveys at mid- and far-IR wavelengths, with the 
aim to parallel optical searches of the deep sky with ob- 
servations at wavelengths where dust is not only far less 
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effective in extinguishing optical light (relevant for esti- 
mating the very uncertain extinction corrections for high 
redshift galaxies, e.g. Meurer et al. 1997), but is also an 
intense source of emission. ISO observations then provided 
an important complementary tool to evaluate the global 
energy output by stellar populations and active nuclei. 

3.1. Overview of the main ISO surveys 

Deep surveys with ISO have been performed in two wide 
mid-IR (LW2: 5-8.5 ^m and LW3: 12-18Aim) and two far- 
IR (A — 90 and 170 /im) wavebands. All surveys are 
performed through repeated raster pointings to achieve 
the best spatial resolution and sensitivity. The diffraction- 
limited spatial resolutions were ^^4.6 arcsec FWHM at 15 
/xm and ~50 arcsec at 100 /im. Mostly because of the bet- 
ter imaging, ISO sensitivity limits in the mid-IR are three 
orders of magnitude deeper in flux density than at long 
wavelengths (0.1 mJy versus 100 mJy). To some extent, 
these different performances are counter-balanced by the 
typical FIR spectra of galaxies and AGNs, which are al- 
most two orders of magnitude brighter at 100 /im than at 
10 /im. We summarize in the following the most relevant 
programs of ISO surveys. 

Five extragalactic surveys with the LW2 and LW3 fil- 
ters have been performed in the ISOCAM Guaranteed 
Time (IGTES), including large-area shallow surveys 
(Si5[Hm] ~ 0.5 — 0.7 mJy) and small-area deep integra- 
tions {Si^llim] ~ 0.1 mJy). A total area of 1.5 square 
degrees have been surveyed in the Lockman Hole and the 
"Marano" southern field, where more than one thousand 
sources have been detected (Elbaz et al. 1999). These two 
areas were selected for their low zodiacal and cirrus emis- 
sions and because of the existence of data at other wave- 
lengths (optical, radio, X). Since the LW2 band at 7/im 
does not sample dust emission in high-z sources and in- 
cludes a large fraction of Galactic stars, we will confine 
our analysis in the following to data in the LW3 15 /im. 

The European Large Area ISO Survey (ELAIS) was 
the largest program in the ISO Open Time (Oliver et al. 
2000a). A total of 12 square degrees have been surveyed at 
15 Atm with ISOCAM and at 90 //m with ISOPHOT (6 and 
1 sq. degrees have been covered at 6.7 and 170 /zm respec- 
tively). To reduce the effects of cosmic variance, ELAIS 
was split into 3 fields of comparable size, 2 in the north 
(Nl, N2), one in the south (SI), plus six smaller areas. 
While data analysis is still in progress, a source list of 
over 1000 (mostly 15 /im) sources is being pubHshed, in- 
cluding starburst galaxies and AGNs (type-1 and type-2), 
typically at z<0.5, with several quasars (including various 
BAL QSOs) found up to the highest z. 

The two ultradeep blank-field exposures by the Hubble 
Space Telescope (one in the North and the other in the 
South, the Hubble Deep Fields, HDF) have promoted a 
substantial effort of multi-wavelength studies aimed at 
characterizing the SEDs of distant and high-z galaxies. 
These areas, including the Flanking Fields for a total of 



~ 50 sq. arcmin, have been observed by ISOCAM at 6.7 
and 15 /im down to a completeness limit of 100 /iJy at 15 
/im. These sensitive ISO surveys have allowed to detect 
dust emission from luminous starburst galaxies up to a 
redshift z = 1.3 (Rowan-Robinson et al. 1997; Aussel et 
al, 1999). In the inner 11 sq. arcmin, the HST provides a 
detailed morphological information for ISO galaxies at any 
redshifts. Thanks to the variety of photometric data and 
an almost complete redshift information available (Aussel 
et al. 1999; Cohen et al. 2000), these surveys are allowing 
the most detailed characterization of the faint IR source 
population. The redshift distributions at the LW3 survey 
limits show an excess number of sources between z=0.5 
and z=1.2, partly an effect of the K-correction as ex- 
plained below. 

Two fields from the Canada-France Redshift Survey 
(CFRS) have been observed with ISOCAM to intermedi- 
ate depths: the '1415-1-52' field (observed at 6.7 and 15 
/im) and the '0302-1-00' field (with only 15 /im data, but 
twice as deep). Studies of ISOCAM sources detected in 
both fields have provided the first tentative interpretation 
of the nature of distant IR galaxies (Flores et al. 1999). 
The LW3 survey displays a redshift distribution similar to 
those of the HDF surveys (see Fig. |l^ below) . 

FIRBACK is a set of deep cosmological surveys car- 
ried out with ISOPHOT, specifically aimed at detecting at 
170 /im the sources of the far-IR background (Puget et al. 
1999). Part of this survey was done in the Marano area and 
in ELAIS Nl, and part in collaboration with the ELAIS 
team in ELAIS N2. This survey is limited by extragalactic 
source confusion in the large ISOPHOT beam (90 arcsec) 
to Si7o > 135 mJy (see for more details Puget et al. 1999 
and Dole et al. 2001). Constraints on the counts below 
the confusion limit obtained from a fluctuation analysis 
of one Marano/FIRBACK field are discussed by Lagache 
& Puget (1999). The roughly 200 sources detected are 
presently targets of follow-up observations, especially us- 
ing deep radio exposures to help reducing the ISO error- 
box and identifying the optical counterparts. Also an ef- 
fort is being made to follow-up these sources with sub-mm 
telescopes (JCMT, IRAM) to derive constraints on their 
redshifts. 

Three lensing galaxy clusters, Abell 2390, Abell 370 
and Abell 2218, have received very long integrations by 
ISOCAM (Altieri et al 1999; Lemonon et al. 1999; Biviano 
et al. 2001). The lensing has been exploited to achieve even 
better sensitivities with respect to ultradeep blank-field 
surveys (e.g. the HDFs), and allowed detection of sources 
between 30 and 100 /iJy at 15 /xm over a total area of 56 
square arcmin (obviously at the expense of an additional 
uncertainty introduced by flux ampliflcation and area dis- 
tortion). The lensing-corrected number counts at 15 /im 
were used by Biviano et al. (2001) to estimate a lower limit 
to the CIRB of 3.3 ± 1.3 nW/m^/sr, close to the upper 
limit by Stanev & Franceschini (1998). 

Ultra-deep surveys in the Lockman Hole and SSA13 
with the LW2 7/im ISOCAM band were performed by 
Taniguchi et al. (1997). The Lockman region was also 
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Fig. 2. Differential counts at Acfr = f5 normalized to the Euclidean law (A^[S'] oc S^'^'^; the differential form is 
preferred here because all data points are statistically independent). The data come from an analysis of the IGTES 
surveys by Elbaz et al. (1999). The dotted line corresponds to the expected counts for a population of non-evolving 
spirals. The short dashed line comes from our modelled population of strongly evolving starburst galaxies, the long- 
dashed one are type-I AGNs. The shaded region at S15 > 10 mJy comes from an extrapolation of the faint 60 /^m 
IRAS counts (Mazzei et al. 2001). 



surveyed with ISOPHOT by the same team: constraints 
on the source counts at 90 and 170 /i are derived by 
Matsuhara et al. (2000) based on a fluctuation analysis. 

3.2. ISO data reduction and analyses 

ISOCAM data needed particular care to remove - in addi- 
tion to the usual photon, readout, flat-field and dark cur- 
rent noises - the effects of glitches induced by the frequent 
impacts of cosmic rays on the detectors (the 960 pixels reg- 
istered on average 4.5 events/sec during the mission). This 
badly conspired with the need to keep them cryogenically 
cooled to reduce the instrumental noise, which implied a 
slow electron reaction time and long-term memory effects. 
For the deep surveys this implied a problem to disentangle 
faint sources from trace signals by cosmic ray impacts. 

To correct for that, tools have been developed by var- 
ious groups for the two main instruments (CAM and 
PHOT), essentially based on identifying patterns in the 
time history of the response of single pixels, which are 
specific to either astrophysical sources (a jump above the 



average background flux when a source falls on the pixel) 
or cosmic ray glitches (transient spikes followed by a slow 
recovery to the nominal background). The most common 
"normal" glitches, induced by cosmic electron impacts and 
lasting only one or two readouts, where the easier to iden- 
tify and remove. Other less frequent impacts by protons 
and alpha particles leave longer-lasting spurious signals, 
from typically several to occasionally one hundreds or 
more readouts. The long integrations adopted for deep 
ISO surveys were needed not only to reduce the instrumen- 
tal noise, but even more to achieve enough redundancy 
(number of elementary integrations per sky position) to 
separate spurious from astrophyical signals in the pixel 
time history. 

A rather performant non-parametric algorithm for 
ISOCAM data reduction is PRETI (Stark et al. 1999), 
exploiting multi-resolution wavelet transforms in the 2-D 
observable plane of the position on the detector vs. time 
sequence). A competely independent parametric method, 
based on a physical model for the detector transients, 
has been devised by Lari et al. (2001), and especially 
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Fig. 3. Integral counts based on the ISOPHOT FIRBACK survey (Sect. 3.1) at AeS = 170 fim (filled circles, from 
Dole et al. 2000) and on the ISOPHOT Serendipitous survey (open squares, Stickel et al. 1998). The shaded region is a 
constraint coming from a fluctuation analysis by Matsuhara et al. (2000). The dashed and dot-dashed lines correspond 
to the non-evolving and the strongly evolving populations as in Fig.^j. The lowest curve is the expected contribution 
of type-I AGNs. The horizontal lines mark the confusion limits (assumed to correspond to 27 beams/source) for 
various telescope sizes: the lines marked "85cm" and "360cm" correspond to the SIRTF and Herschel/FIRST limits 
for faint source detection, assumed diffraction-limited spatial sampling. The confusion limit scales up on the y-axis 
proportionally to the square of the diameter of the primary reflector. 



taylored for the shallow ELAIS integrations. Independent 
tools have been developed by Desert et al. (1999) and by 
Oliver et al. (2000a). 

The PRETI and LARI detection algorithms have been 
tested by means of Monte Carlo simulations including 
all artifacts introduced by the analyses. Such simulations 
have been performed on real datasets, including both a 
long staring observation of more than 500 readouts (Elbaz 
et al. 1998) and the deep survey frame itself (Lari et al. 
2000). Test sources with known fluxes were introduced 
with their PSF and model transients, against which the 
detection algorithm has been tested. With these simula- 
tions it is has been possible to control as a function of the 
flux threshold: the detection reliability, the completeness, 
the Eddington bias and photometric accuracy (~10-20% 
when enough redundancy was available, as in the CAM 
HDFs and IGTES ultradeep surveys). The PRETI and 
LARI methods have been applied in particular to the HDF 
North dataset (Aussel et al. 1999, C. Lari, private com- 
munication), and they showed excellent agreement down 
to the faintest fluxes. 

The astrometric accuracy is of order of 2 arcsec for 
deep highly-redundant and properly registered images, al- 
lowing relatively easy identification of the sources (Aussel 



et al. 1999, 2001; Elbaz et al. 2001; Fadda et al. 2001). 
For example, among the complete ISOCAM sample of 41 
galaxy identifications in the HDF North studied by Aussel 
et al. (1999, 2001), all have an optical counterpart brighter 
than 1=23 within 3", and only one source appears to be 
confused by the presence of more than one optical galaxy 
in the ISO errorbox. The quality and reproducibility of 
the results for the CAM surveys is also proven by the 
good consistency of the counts from independent surveys 
(see Fig.[H below). 

A reliable reduction of the longer wavelength 
ISOPHOT observations proved to be more difficult. The 
170^*™ counts from PHOT C200 surveys are 90% reli- 
able above the 5-cr confusion limit S170 ~ 130 mJy, and 
required relatively standard procedures for baseline cor- 
rections and " de-glitching" . Quite more severe are the 
noise problems for the ClOO channel (60 to 90/iTO, which 
would otherwise benefit by a better spatial resolution), 
preventing so far to achieve significantly better sensitivi- 
ties than IRAS. The ClOO PHOT survey dataset is still 
presently under analysis (C. Lari and G. Rodighiero, work 
in progress). 
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Fig. 4. Differential counts from the ISOPHOT FIRBACK 
survey at Acff — 170 /im, from Dole ct al. (2001), compared 
with the prediction of the reference model. As in Fig. ||, the 
dotted line corresponds to the non-evolving spiral popula- 
tion, the short-dashed to evolving starbursts, long-dashes 
are type-I AGNs. 

3.3. Mid-IR and far-IR source counts from ISO surveys 

IR-selected galaxies have typically red colors, partly be- 
cause of extinction by dust responsible for the excess IR 
flux. When found at substantial redshifts, these sources 
are also quite faint in the optical. For this reason the 
redshift information is presently available only for limited 
subsamples (e.g. the HDF North and CFRS samples). In 
this situation, the source number counts provide crucial 
constraints on the evolution properties. 

Particularly relevant information comes from the mid- 
IR surveys based on the ISOCAM 15 urn LW3 filter, be- 
cause they include the faintest, most distant and most nu- 
merous ISO-detected sources with reliable identifications. 
To cover with LW3 a wide dynamic range in flux with good 
source statistics, Elbaz et al. (1999) performed a variety 
of surveys with sky coverages decreasing as a function of 
the flux limits. Including ELAIS and the IRAS data, the 
range in fluxes reaches four orders of magnitude. 

The differential counts (normalized to the Euclidean 
law dN oc S~''^-^dS) based on data from seven independent 
sky areas, shown in Fig. ^, reveal a remarkable agreement. 
Of the various samples considered by Elbaz et al., only 
sources in flux bins for which the survey was better than 
80% complete were used (for a total of 614 sources). 

At fluxes fainter than 1 mJy the contamination by 
stars is of the order of only few %, while at brighter fluxes 
in the Lockman Shallow Survey it reaches ~ 10%, and 
further increases at increasing flux. 

In addition to the data reported by Elbaz et al., the 
shaded region at S15 > 5 mJy in Fig. ^ corresponds to an 



estimate of the extragalactic counts by Mazzei et al. (2001) 
overriding the problem to account for the large fraction 
of galactic stars at these bright fluxes. This estimate is 
based on the 60 /im IRAS galaxy counts translated to the 
LW3 band by using ISOCAM photometry of a complete 
sample of faint IRAS sources in the North Ecliptic Pole 
(Aussel et al. 2000). This evaluation of the bright 15 fim 
counts helps to constrain the level of non-evolving galaxies 
and the normalization of the local luminosity function (see 
Sect. 5.1), untill a systematic identification of the ELAIS 
catalogues will be available. 

The combined 15 fim differential counts display vari- 
ous remarkable features (Elbaz et al. 1999): a roughly eu- 
clidean slope from the brightest fiuxes sampled by IRAS 
down to S15 ~ 10 mJy; a fast upturn at S15 < 3 mJy, 
where the counts increase as dN oc S~^'^dS to S15 ^ 0.4 
mJy; and finally the evidence for a flattening below S15 ^ 
0.3 mJy (where the slope becomes quickly sub-Euclidean, 
dN oc S~^dS). Note that the sudden change in slope and 
the faint flux convergence is supported by 3 independent 
surveys. 

The areal density of ISOCAM 15 /im sources at the 
limit of ~50/iJy is ~5 arcmin^^. If we consider that 
the diffraction-limited diameter of a point-source is ~ 50 
arcsec^ and for a slope of the counts (3 ^ —2, this den- 
sity is close to the ISO confusion limit at 15 /im of ~ 0.1 
sources/areal resolution element, or 7 / arcmi'n? in our case 
(Franceschini 2000). Confusion will likely remain a limita- 
tion for the NASA's SIRTF mission, in spite of the mod- 
erately larger primary collector. 

Obviously, far-IR selected samples are even more se- 
riously affected by confusion. The datapoints on the 
170/<m integral counts reported in Fig. § come from the 
FIRBACK survey (Dole et al. 2001), while Fig. | shows 
the same counts in differential units. Similarly deep ob- 
servations at 90, 150 and 180 /im are reported by Juvela, 
Mattila & Lemkc (2000). These surveys have essentially 
attained at the ISOPHOT confusion limit. Additional con- 
straints at slightly fainter fluxes have been attempted us- 
ing background fluctuation analyses (Lagache & Puget 
1999; Matsuhara et al. 2000, see shaded region in Fig. 

I). 

Shorter wavelength ISOPHOT ClOO observations 
could in principle benefit by a less severe confusion lim- 
itation. Preliminary results of ISOPHOT ELAIS surveys 
at Aoff — 90 /im (Efstathiou et al. 2000), as well as counts 
derived from the IRAS 100 /im survey, are reported in Fig. 
^, showing nice agreement in the overlap fiux range. 

4. GALAXY SOURCE COUNTS FROM 
SURVEYS AT MILLIMETRIC 
WAVELENGTHS 

Surveys in the sub-millimeter offer a unique advantage 
to naturally generate volume-limited samples from flux- 
limited observations. This property is due to the peculiar 
shape of galaxy spectra, with an extremely steep slope 
from 1 mm to 200 /im [roughly L{v) oc z^'^'^, Andreani & 
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Fig. 5. Integral counts at Acff = 90 /im based on the IRAS 
(open circles) and the ELAIS 90 /im (filled circles) surveys, 
as reported by Efstathiou et al. (2000). Line symbols as in 
Fig. ^. The lines marked "60cm" and "360cm" correspond 
to the ISO and Herschel/FIRST limits for faint source 
detection. 

Franceschini (1996)]. Then, as we observe in the sub-mm 
galaxies at larger and larger redshifts, the selection wave- 
band in the source rest-frame moves to higher frequencies 
along a steeply increasing spectrum, and the correspond- 
ing K-correction almost completely counter-balances the 
cosmic dimming of the flux, for sources at z > 1 and 
up to z ~ 10 (Blain & Longair 1993). A further related 
advantage of sub-mm surveys is that local galaxies emit 
very modestly at these wavelengths. Altogether, a sen- 
sitive sub-mm survey will avoid local objects (stars and 
nearby galaxies) and will select preferentially sources at 
high and very high redshifts: a kind of direct picture of 
the high-redshift universe, impossible to obtain at other 
frequencies. 

Important discoveries have come from the operation 
of a powerful array of bolometers (SCUBA) on JCMT, 
due to a combined effect of sensitivity, large multiplexing 
capability, an efficient 15m primary mirror, and an op- 
eration site allowing to observe at short mm wavelengths 
(A = 850 [mi). In such way, SCUBA on JCMT has allowed 
to resolve more than 20% of the long-A CIRB background 
into a population of faint distant, mostly high-z, sources. 

Various groups have used SCUBA for deep cosmolog- 
ical surveys. Smail et al. (1997, 1999) have exploited 7 
distant galaxy clusters as cosmic lenses, obtaining a sam- 
ple of 17 (3(t) sources with Sgso > 6 mJy. Hughes et al. 
(1998) published a single very deep image of the HDF 
North containing 5 (4cr) sources at Sgso > 2 mJy. Lilly 
et al. (1999) and Bales et al. (2000) have published sam- 
ples including '^30 (Str) sources to 3.5 mJy. Scott, Dunlop 
et al. (2001, in prep) have recently completed a survey 



in the Lockman and ELAIS N2 regions and detected ap- 
proximately 20 (4cr) sources. In addition to deep sub-mm 
surveys (Barger et al. 1998, 1999), Barger, Cowie and co- 
workers have carried out an extensive program of follow-up 
of SCUBA sources with optical telescopes. Other system- 
atic identification efforts have been attempted by Ivison 
et al. (2000, see also Sect. 6.3). 
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Fig. 7. Integral counts at Aofr = 1200 /im by Bertoldi et 
al. (2001). See also caption to Fig. ^. 

Each of these deep integrations required many tens 
of hours each of especially good weather, which meant 
approximately 20% of the JCMT observatory time since 
1997. In spite of this effort, the surveyed areas (few tens 
of arcmiv?) and number of detected sources (from 100 to 
200 sources) are modest, which illustrates the difficulty 
to work from ground at these wavelengths. Much fewer 
sources have been detected in the 450 yuva channel, for 
which the atmospheric transmission at JCMT is usually 
poor. 

The extragalactic source counts at 850 /im, reported in 
Fig. ^, show a dramatic departure from the Euclidean law 
\dN cx dS in the crucial flux-density interval from 2 
to 10 mJy], a clear signature of the strong evolution and 
high redshift of SCUBA-selected sources. 

A new bolometer array (MAMBO) has become re- 
cently operative on the largest existing mm telescope, the 
IRAM 30m. A large survey at XeS — 1-2 mm of 3 fields 
with a total area of over 200 arcmin^ to a flux limit of 
few mJy has been performed with MAMBO (Bertoldi et 
al. 2000): preliminary galaxy counts and evaluations of 
the redshift distributions are reported by Bertoldi et al. 
(2001), both confirming the SCUBA results at 850 /im 
(see Fig. 0). 
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Fig. 6. Integral counts at AeS = 850 /um. The filled squares are from Blain et al. (1999), the open squares from Barger, 
Cowie and Sanders (1999). See caption to Fig. || for meaning of the lines. 



5. INTERPRETATIONS OF FAINT IR/MM 
GALAXY COUNTS AND 
REDSHIFT-DISTRIBUTIONS 

5.1. Predictions for non-evolving source populations in 
the mid-IR 

A zero-th order interpretation of data on deep counts is ob- 
tained by comparing them with the expectations of models 
assuming no-evolution for cosmic sources. When referred, 
in particular, to the mid-IR galaxy counts of Fig. ^, this 
calculation has to account in detail for the effects of the 
complex mid-IR spectrum of galaxies (including strong 
PAH emission and silicate absorption features, see Fig. ^ 
and |l^ below) in the K-correction factor relating observed 
flux S and rest-frame luminosity L: 



L,K{L,z) 
4ndl ' 



(5) 



where c^l is the luminosity distance Q and K{L, z) = 
(i+z)L[iy(i+z)] K-correction. The differential number 
counts (sources/unit flux interval/unit solid angle) at a 
given flux S write as: 



dN 

Is 



dz — — 

dz 



dV d\ogL{S,z) 



dS 



p[L{S,z),z] 



(6) 



The luminosity distance can be evaluated for SIa = and 
any value of fim from the usual Mattig relation di, — z[l + 
z{l-nm/2)/{Vl + nmZ + l + fln,z/2)]. Foi + = 1, then 
rfL = (1 + z) J^[{1 + z')'\l + n^z') + z'{2 + z')nxr^^^dz. 



where p[L{S,z),z] is the epoch-dependent luminosity 

function and dV/dz = ^prj^^^j^^ is the differential vol- 
ume element. The moments of various order of the distri- 
bution dN/dS of eq.(ph provide respectively the integral 
counts A^(> S) — J ^dS, the contribution to the back- 
ground intensity — J dSi,, and to the back- 
ground fluctuations < [61{iy)]'^ / dS^. 

Taking into account the system transmission function 
T(A), the K-correction is more appropriately written as: 



K{L,z)^ 



jl + z) /,^;dA if) r(A) mi + z)] 

/,';dA {^) TiX)L[u,z = 0] 



(7) 



Fig. H illustrates the effects of the mid-IR emis- 
sion/absorption features on K-correction, for different 
mid-IR spectra (an inactive spiral [dotted line] , a M82-like 
starburst [continuous line], both observed with ISOCAM 
LW3) and different filters (the starburst M82 observed 
with LW2 [dashed fine]). The detailed spectral response 
functions of the various filters, as given in the ISOCAM 
Observer Manual (ESA Publications), have been used in 
eq. (^. The effect of the prominent emission features on 
the counts is particularly important in the wide LW3 (12- 
18 /im) filter, as they imply an enhanced sensitivity to 
sources at 0.5 < z < 1.3, and a diminished one out- 
side this redshift interval. Correspondingly, the LW3 se- 
lection provides improved capability to study in detail the 
evolution of galaxies between z — 0.5 and 1, a redshift- 
interval that current cosmogonic models (e.g. Kauffmann 
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Fig. 8. K-correction as a function of the source redshift 
(eqj^) for different filters and source spectra. Dotted line: 
inactive spiral spectrum observed with the LW3 filter; con- 
tinuous line: M82-like spectrum with LW3; dashed line: 
M82 spectrum observed with LW2. Note that the spiral 
spectrum (dotted line) implies the strongest K-correction 
at z ~ 1 because of the lack of hot-dust emission depress- 
ing the rest-frame LW3 flux compared with the redshifted 
PAH-dominated emission. 



& Chariot 1998) predict to be critical for the formation of 
structures. 

The second ingredient for the no-evolution prediction 
is the local galaxy luminosity function (LLF). In the mid- 
IR, LLPs have been published by Rush et al. (1993), Xu 
et al. (1998) and Pang et al. (1998), all based on the 12 
fim all-sky IRAS survey. Unfortunately, in spite of the 
proximity of the CAM LW3 and IRAS 12 ^m bands, our 
knowledge of the 15 /im LLP is still somehow uncertain, 
because of: a) uncertainties in the IRAS 12 /im photome- 
try; b) the effect of local inhomogeneities, particularly the 
local Virgo super-cluster, in the shallow IRAS survey; and 
c) the uncertain flux conversion between the IRAS and 
CAM-LW3 bands (Alexander & Aussel 2000). Because of 
(a) and (b), the Rush et al. LLP determination was af- 
fected by an improper flux normalization and a too steep 
faint-end slope (e.g. inconsistent with the IRAS 60 /im 
LLP). These various effects have been discussed by Pang 
et al. (1998) and Xu et al. (1998): a re-calibrated 12 /im 
luminosity function based on these analyses is reported in 
Pig. ^. Open squares at L12 > 10^" L0 come from Pang 
et al. (1998). Por lower luminosity values the Fang et al. 
LLP persisted to show excess density, and for this reason 
we used here the determination by Xu et al. (1998, we 
neglect the small correction from 15 to 12 /im, given the 



flat shape of the LLF here) . The shallower low- luminosity 
slope of this new LLF determination is in particular con- 
sistent with the 60 /im one (small filled squares in Fig.^ 
see also Sect. 5.3 below). We fit these data (dotted line in 
Pig. 1^) with an analytic form 

p(L) = /,, xL(i--) x(l + i/L,/6)-^ 

similar to the one suggested by Rush et al., but with dif- 
ferent slopes (a=1.15 vs. 1.7, b=3.1 vs. 3.6) and different 
normalizations. 

The dotted line in Fig. ^ corresponds to the pre- 
dicted counts assuming a non-evolving population using 
this best-estimate LLP. The correction to the CAM LW3 
band is made by adopting a 12 /im to LW3 flux ratio which 
is a function of the 12 /xm luminosity: for the less luminous 
objects the ratio is based on the observed mid-IR spec- 
trum of quiescent spirals, while for the highest luminosity 
galaxies the ratio is the one expected for ultraluminous IR 
galaxies, and for intermediate objects is close to a typical 
starburst spectrum like the one of M82 (continuous line 
in Fig. |l^). The 15 to 12 /xm flux ratio is then assumed 
to increase continuously with luminosity, the 15 /im flux 
being increasingly dominated by the starburst emission at 
increasing L (see Sect. 5.3 for details). Note that in the 
absence of evolution, different values for the cosmological 
parameters have negligible influence on the no-evolution 
expectation. 

It is clear that the no-evolution prediction, even tak- 
ing into account the effects of the PAH features on the 
K-corrections, falls quite short of the observed counts at 
fluxes fainter than a few mJy. Also the observed slope 
in the 0.4 to 4 mJy flux range {dN[S] cx S~^^°-^dS in 
differential units) is very significantly different from the 
no-evolution prediction {dN{S) cx S^'^dS). The extrapo- 
lation to the bright fluxes is instead consistent, within the 
uncertainties, with the IRAS 12 /im counts with a slope 
close to Eclidean. 

Another clear sign of a serious mis-match is provided 
by the comparison of the no-evolution prediction with the 



observed redshift distributions in Fig. 10, where the model 
keeps a large factor below the observed peak at 2; 0.85. 

5.2. Evidence for a strongly evolving population of 
mid-IR galaxies: the reference model 

A first robust conclusion of the previous Section was 
the need for an evolving source population dominating 
the counts. The shape of the differential counts shown 
in Pig. 1^ contains some indications about the properties 
of the evolving populations. In particular the almost flat 
(Euclidean) normalized counts extending from the bright 
IRAS fluxes down to a few mJy, followed by the sudden 
upturn below, suggests that is not likely the whole popu- 
lation of IR galaxies that evolve: in this case and for the 
adopted IR galaxy LLP, the super-Euclidean increase in 
the counts would appear at brighter fluxes and not be as 
abrupt as observed. The observed behaviour is more con- 
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Fig. 9. Galaxy LLF's at 12 /ini from Fang et al. (1998, red open squares) and adapted from Xu et al. (1998) in 
the low-luminosity regime. A comparison is made with the IRAS 60 fim LLF by Saunders et al. (1990, small filled 
squares). Black circles are an estimate, based on the Rush et al. (1993) catalogue, of the contribution to the 12 /im LLF 
by active galaxies (including type-I AGNs [red long-dashed line] and type-II AGNs plus starbursts [green short-dash 
line]). Type-II AGNs and starbursts are included in the same population on the assumption that in both classes the IR 
spectrum is dominated by starburst emission. The dotted line is the separate contribution of normal spirals, while the 
continous line is the total LLF. In our reference model (Sect. 5.2), only type-I AGNs and active galaxies dominated 
by starburst emission (short-dashes) are allowed to evolve. 



sistent with a locally small fraction of IR galaxies to evolve 
with high rates back in cosmic time. 

We have reproduced the IR counts in Fig. || with the 
contribution of three population components character- 
ized by different physics and evolutionary properties. The 
main contributions come from non-evolving normal spi- 
rals (with a 12 /im LLF as the dotted line in Fig. ||), and 
from a fast evolving population (LLF as the green dashed 
line in Fig. ^). The evolving population includes type- 
II AGNs and starburst galaxies, with the idea that for 
both classes the IR spectrum may be dominated by star- 
burst emission. A third component considered are type- 
I AGNs, whose LLF is the long-dashed line in Fig. ^. 
Based on optical and X-ray surveys (e.g. Franceschini 
et al. 1994b), the latter are assumed to evolve in lumi- 
nosity as L{z) — L(0) X (1 -I- z)^ up to z = 1.5 and 
L{z) — L{0) x2.5'^ above. This different treatment of type- 
I and II AGNs is not necessarily in contradiction with the 
Unified AGN Model: for type-II objects the inclined dusty 



torus implies a self-absorbed emission in the mid-IR by the 
AGN, and an overall IR spectrum dominated by a circum- 
nuclear starburst. 

A similar multi-population modelling was proposed 
time ago by Danese et al. (1987) to explain the IRAS and 
faint radio counts, and more recently by Roche & Bales 
(1999) and Xu et al. (2001). 

As shown in Fig. ^, the fraction of the evolving star- 
burst population is assumed to be ~ 10 percent of the 
total, roughly consistent with the local observed fraction 
of interacting galaxies. The quick upturn in the counts 
then requires a strong increase with redshift of the av- 
erage emissivity of the evolving population to match the 
peak in the normalized counts around S15 ~ 0.5 mJy. 

In the presence of strong evolution, the fit to the ob- 
served counts has a sensible dependence on the assumed 
values for the geometrical parameters of the universe. For 
a zero-A open universe (ilm = 0.3, J^a — 0), a physi- 
cally plausible solution would require a redshift increase 
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Fig. 10. Redshift distribution of 15 /im sources with S15 > 
120 /iJy in the HDF North (continuous red histogram; 
Cohen et al. 2000; Aussel et al. 2001), compared with our 
best-fit evolutionary model (continuous line) . The dashed 
green histogram is the z-distribution of sources in the 
CFRS field reported by Flores et al. (1999). The dotted 
line at the bottom corresponds to the no-evolution predic- 
tion. 

of the comoving density of the starburst sub-population 
as p{L[z\, z) = po{Lq) X (1 + z)^ and of their luminosity as 
L(z) = Lq X (1 + zY for z < Zbreak, with p and L constant 
above, and Zbrcak = 0.9. These are quite extreme evolu- 
tion rates, if compared for example with those observed 
in optically-selected samples of merging and interacting 
galaxies (e.g. Le Fevre et al. 2000). 

The inclusion of a non-zero cosmological constant, with 
the consequent increase of the cosmic timescale and vol- 
umes at z ~ 1, tends to make the best- fitting evolution 
rates less extreme. For flai = 0.3, f^A — 0.7, a best-fit to 
the counts requires: 

p{Li2[z],z) = pa{Li2) X (1 + z)'' Z < Zbroak 
p{Li2[z],z) = Po{Li2) X (1 + Zbrcak)** Zbrcak < Z < Z,„ax 
Li2(z) = L12 X (1 + z)^-^ Z < Zbroak 

Ll2(z) = L12 X {1 + Zbrcak)^'* Zbroak < Z < Zmax (8) 
with Zbroak = 0.8 and Zmax = 3.7. 

Note that, although there are margins for variations 
of the relative weight for the number density and lumi- 
nosity evolution, only a combination of the two provides 
credible solutions. Assuming for example evolution only 
in source number density or only in L would require un- 
plausibly high evolution rates for the evolving population 



(p oc [1 -I- z]^^, or alternatively L[z] oc [1 + z]^). Such ex- 
treme solutions would also encounter problems when fit- 
ting the radio or far-IR counts as counterpart to the 15 
/xm population counts. 

The vast majority (> 90%) of the sources in the HDFN 
and CFRS 1415 ISO surveys have spectroscopic redshifts, 
and for the remaining objects photometric redshifts are 
easy to estimate. The redshift distributions D{z) for the 
HDF North (Cohen et al. 2000; Elbaz et al. 2001; Aussel 
et al. 2001) and CFRS 1415 surveys (Flores et al. 1999) 
are compared in Fig. ^ with our best-fitting model, which 
appears to properly account for these data. They set strin- 
gent limits on the rate of cosmological evolution for IR 
galaxies above z ~ 1, and force it to level off to avoid 
exceeding the observed D{z) on the high-z tail. Also, a 
consequence of the fast evolution observed at z < 1 is 
that the observed CIRB intensity is quickly saturated at 
moderate z, and requires again the evolution rate to turn 
over at z > 1. 

Unfortunately, apart from these constraints on the 
high-z IR emissivity, ISO surveys do not offer an accu- 
rate sampling of the hidden SF in the interval 1 < z < 2, 
which will only be possible with the longer- wavelength sur- 
veys by SIRTF at XeS = 24 /xm and by the Herschel Space 
Observatory (formerly FIRST) in the far-IR, probing dust 
and PAH emission at z > 1. 

In our present evolutionary scheme, any single galaxy 
would be expected to spend most of its life in the quies- 
cent (non-evolving) phase, while being occasionally put by 
interactions in a short-lived (few to sevral 10*" yrs) star- 
bursting mode. The cosmological evolution characterizing 
this phase may simply be due to an increased probabil- 
ity in the past to find a galaxy during such an excited 
mode. The density evolution in eq. (^) scales with red- 
shift approximately as the rate of interactions due to a 
simple geometric effect following the increased source vol- 
ume density. The luminosity evolution may be interpreted 
as an effect of the larger gas mass available to form stars 
at higher z. 

5.3. A panchromatic view of IR galaxy evolution 

Deep surveys at various IR/sub-mm wavelengths can be 
exploited to simultaneously constrain the evolution prop- 
erties and broad-band spectra of faint IR sources. We re- 
port in this Section a comparison of the 15 pm survey 
data with those coming from longer-wavelength surveys, 
in particular the IRAS 60 pm, the FIRBACK 90 and 170 
/im, and the SCUBA 850 pm data, which are the deep- 
est, most reliable available at the moment. Information on 
both number counts and the source redshift distributions, 
whenever available, were used for these comparisons. 

Further essential constraints, providing the local 
boundary conditions on the evolutionary histories, are 
given by the multi-wavelength local luminosity functions. 
In addition to the 12 and 15 pm LLF's discussed in Section 
5.1, the galaxy LLF is particularly well known at 60 pm 
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Fig. 11. Local luminosity function of galaxies at Aofr — 
60 (filled squares) by Saunders et al. (1990). The lines 
are the best-fit model predictions (upper continuous line: 
total; dotted: quiescent population; lower continuous: star- 
bursts). 

after the IRAS all-sky survey and the extensive spectro- 
scopic follow-up (Saunders et al. 1990), and is illustrated 
in Fig. ^ Dunne et al. (2000) also attempted to constrain 
the galaxy LLF in the millimeter based on mm observa- 
tions of complete samples of IRAS 60 /im galaxies. The 
results are shown in Fig. 

As previously mentioned, the properties of LLF's ob- 
served at various IR/sub-mm wavelengths can be ex- 
plained by assuming that the galaxy IR SED's depend on 
bolometric luminosity. The comparisons of LLF's made in 
Figs. |9| and |lj show that the 60/im LLF has a much flat- 
ter power-law shape at high-luminosities compared with 
both the mid-IR and millimetric LLF's. This is clearly an 
effect of the spectra for luminous active galaxies showing 
excess 60 /xm emission compared with inactive galaxies, as 
also illustrated by the luminosity-dependence of the IRAS 
far-IR colours. We defer to Franceschini (2000, Sect. 6.6) 
for further discussion on this effect. 

Consequently, we have modelled 

the redshift-dependent multi-wavelength LLF's of galax- 
ies by assuming spectral energy distributions dependent 
on luminosity, with spectra ranging from those typical of 
low-luminosity inactive objects, to those peaked at 80 /im 
of luminous and ultra-luminous IR galaxies as previously 
described. 

We have taken as reference for our multi-wavelength 
LF the one by IRAS at 12 ^m (p(Li2)) discussed in Sect. 
5.1. This is transformed to longer wavelengths according 
to spectral energy distributions which vary according to 



Fig. 12. Local luminosity function of galaxies at Aoff = 
850$mum (open squares) by Dunne et al. (2000), com- 
pared with the IRAS 60 fim LLF by Saunders et al. (1990). 
Note the completely different slopes of the two. 

the value of L12. The assumption was that for L12 < Lmin 
the spectrum is that of an inactive spiral (Siowl^^], the 
lower dotted line in Fig. while for L12 > Lmax it is a 
typical ULIRG spectrum (Shigh [i^] , top line in Fig. ^ . For 
intermediate luminosity objects, the assumed SED S(j^) is 
a linear interpolations between the two: 

= (9) 

Slo^{iy){l0gLi2 - logLmin) + 5'higli (t^) (log L12 - logLmjn) 

logL 

max 

- log iinin 

From <S'(i^)|Li2 it is immediate to compute the luminos- 
ity at any frequencies and the K-correction from eq. (^, 
taking into account the detailed filter response functions. 
The multi- wavelength LLFs at any wavelengths A are eas- 
ily computed from the relation 

P(L.|L,2,.)^P(L,2,.)(^)". (10) 

Good fits to the multi- wavelength LF's are found by set- 
ting Lmin = 2 10^ Lq and L^ax = 10^^ Lq. 

For the evolving active starburst galaxies we adopted 
both a single average spectral energy distribution (inde- 
pendent on luminosity) and a luminosity-dependent spec- 
tral shape as discussed for the non-evolving population. 
For simplicity and for a better controlled parametrization, 
our best-fit model for the active starburst population as- 
sumes the single spectrum solution. 

If we adopt as representative for the starburst spec- 
trum the IR SED of the ultra-luminous galaxy Arp 220, 
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would predict that of order of 4% of the FSS could be 
starbursts at z > 0.4. 

It is evident that, with various degrees of significance 
which depend on the survey depths, all the observed long- 
wavelength counts require a substantial increase of the IR 
volume emissivity of galaxies with redshift. 




A(/iim) 



Fig. 13. Our adopted range of mid- to far-IR spectra 
of galaxies. The lower dotted line corresponds to a low- 
luminosity inactive spiral (Siow[H in eq. |^), the top dotted 
line is typical of ULIRGs (Shigh[H)- The intermediate one 
(cyan continuous line) corresponds to our adopted spec- 
trum for the starburst population: this spectrum is similar 
to the one of the prototype star-forming galaxy M82 (in 
the range from 5 to 18 /im it is precisely the ISO CAM 
CVF spectrum of M82). 



the consequence would be that all far-IR counts (and the 
CIRB intensity, see Elbaz et al. 2001) would be exceeded 
by substantial factors. On the contrary, if we assume for 
the IR evolving sources a more typical starburst spectrum 
(see continuous line in Fig. |l^, which is similar to those 
of M82 and other luminous starbursts observed by ISO), 
then most of the observed properties of far-IR galaxy 
samples (number counts, redshift distributions, luminosity 
functions) are appropriately reproduced. 

Best-fits to the counts based on our reference model 
and adopting a spectral template as in Fig. |l3|, are given 
in Figs. I, |, |, I and |. Fig. |l| reports the fit to the 60 
/xm differential counts derived from various IRAS surveys, 
while Fig. |l^ compares our model prediction with the z- 
distribution for faint galaxies selected from the IRAS Faint 
Source Survey at 60 /zm by Oliver et al. (1996, see also 
Saunders et al. 2000), with a flux limit of Seo = 200 mJy. 
In addition to the 60 /im counts, the model clearly repro- 
duces the observed z-distribution up to z = 0.4, while it 
predicts a somewhat excess fraction of higher-z sources. 
We do not interprete this as to necessarily be a prob- 
lem for the model, since due to the large IRAS errorbox, 
the most distant and optically faint sources could have 
been systematically mis-identified with brighter more lo- 
cal galaxies falling by chance in the errorbox: our model 



Fig. 14. Differential counts at 60 fim from IRAS surveys, 
normalized to the Euclidean law 600 S~'^'^ {sr~^ J~^), ver- 
sus the 60 /im flux in Jy. Meaning of the lines as in previous 
figures. 




Fig. 15. Predicted redshift distribution for the 60 /im 
IRAS Faint Source Survey identified by Oliver et al.(1996). 
The survey depth and area are indicated. Meaning of the 
lines as in the previous figures. 

We report in Fig. |l^ the evolutionary 15 /xm luminos- 
ity function of IR sources in the critical redshift interval 
from z=0 to 1. The evolution of the high- luminosity end is 
partly driven by the type-I AGN population, which dom- 
inates the LF above L15 = 10^^ Lq. Note that p{L,z) is 
required to exceed at z ~ 1 the local galaxy density (as 
evident in particular in the luminosity interval L15 ~ 10*" 
to L15 ~ lO"'^" Lq), which we interprete as an effect of 
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merging, implying an increase of the comoving number of 
objects in the past. 

The evolution pattern of the global LF is clearly lumi- 
nosity dependent: the largest increase happens at L15 ~ 
10^^ Lq, whereas at very low and very large luminosities 
the evolution is lower. The model predicts that the evolu- 
tion of galaxies with L15 > 10^^'^ Lq (the Ultra-Luminous 
Galaxies, ULIRGs) is lower than that of the Luminous IR 
galaxies (LIRGs) around 10^^ Lq. The assumption of sim- 
ilar evolution for LIRGs and ULIRGs would tend to imply 
exceeding the observed fraction of z > 0.4 galaxies in the 
60 fiui IRAS Faint Source Surveys and the North Ecliptic 
Pole survey (Ashby et al.l996, Aussel et al. 2000). 




10' 10' 10" 10i» 10" 10'2 

log Li=,„ (LJ 



Fig. 16. Evolution of the 15 /im luminosity function of IR 
sources, following our best-fitting reference model. The 
lower continuous and dotted lines are the contribution of 
AGN and of normal galaxies to the LLF. The three con- 
tinuous lines labelled z=0, 0.4 and 1 are the global (AGN 
-|- galaxy) LFs at the corresponding redshifts. 



5.4. Alternative evolution patterns 

The available data at long— A on faint galaxies, although 
constraining the evolution pattern, do not allow to deter- 
mine it univocally. We have attempted, in particular, evo- 
lutionary schemes considering a single evolving population 
with combined luminosity/density evolution, as an alter- 
native to the previously discussed two-population model. 



5.4.1. Luminosity-independent evolution for a single 
population 

This model allows the LF for the whole population of IR 
galaxies to evolve in both luminosity and volume density 
following a pattern as in eq. (^, with reduced evolution 
rates: p(Li2, z) oc {1 + z)'^-^, L12 oc (1 + z)^, Zbroak = 0.8, 
Zmax = 3.5, and Lmin = 3 10^, Lmax = 10^^ in eq. (pi), 
Siow(j^) and Shigh('^) similar to those plotted in Fig. [13| . 

Forcing it to best-fit to the 15 /im counts and D(z), the 
model appears to be too rigid when fitting data at longer 
wavelengths. In particular, the predicted 60 /um counts 
show a fast convergence already below Seo ~ 0.1 Jy and a 
very modest contribution to the CIRB, while at the same 
time the 170 /im counts would be very steep and exceeding 
the observations. The predicted 850 /im counts are also 
steeper than observed. 

5.4.2. Luminosity-dependent evolution rates 

Under this scheme, the evolution rates are assumed to 
depend linearly on luminosity: p{Li2,z) oc (1 -f zy^^^^\ 
with j{Li2) = J X log{Li2/Li)/log{L2/Li), and L12 oc 
(1 + z)'=(^i^), with fc(Li2) =Kx log{Li2/Li)/log{L2/Li) 
[with J, K, Li and L2 as adjustable parameters]. This 
law is intended to add freedom to the single-population 
model, by considering that galaxies with higher L12 (hence 
more "active") may evolve with faster rates with respect 
to lower-Li2 (less "active") objects. 

As a matter of fact, the larger freedom for this model 
is of no help to improve the quality of the fits (the param- 
eters Li and L2 tend to assume low values, i.e. to go back 
to the luminosity-independent solution). 

Altogether, evolutionary schemes alternative to our 
reference two-population model appear to provide worse 
fits to the multi- wavelength combined IR data. Our previ- 
ously discussed evolution model provides the best statisti- 
cal fit of the data with a plausible physical interpretation 
(but certainly not the only acceptable solution). 

5.4.3. Variations in the model population combination 

Finer variations of the relative fractions of sources belong- 
ing to the different physical components (AGNs, normal 
spirals and evolving active galaxies), with respect to the 
values adopted in our best-fit model, are clearly allowed 
by the present data. However, the 15 /Ltm number counts, 
local LF and redshift distributions, in particular, together 
with looser constraints from longer-wavelength observa- 
tions, allow rather narrow margins to such variations. For 
example, if the normalization of the evolving population 
is incresed and that of the non-evolving one decreased to 
still keep fitting the 15 fim local LF, we may obtain steeper 
counts in Fig. ^ in the flux interval from 10 to 0.5 mJy 
(with a suitable change of the evolution rates). This, how- 
ever, would start forcing the fit to the observed counts in 
this flux range, and would also tend to spoil those to the 
longer- A data (Figs. 3, 11, 14). 
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Note also that our modelling of the AGN component 
cannot be but rough at this stage. Our adopted proce- 
dure treats type-I and type-II AGNs separately. The lat- 
ter are simply included in a single evolving population 
together with the active starbursts. This reflects our view 
that dust-extinguished AGN and starburst emissions hap- 
pen concomitantly in the same sources during the "ac- 
tive" phase and reflects our present inability to quantita- 
tively disentangle the two. This concomitant "activity" is 
now revealed by a variety of facts, including the hard X- 
ray observations of ISO sources reported by Fadda et al. 
(2001, see also Elbaz et al. 2001), and some evidence that 
the brightest sub-mm objects are associated with AGNs 
(Ivison et al. 2000a). 

On the contrary, our modelling of the type-I AGNs is 
simple and robust, and exploits the 12 /im LLF by Rush et 
al. (1993) and the evolution rate found from optical and 
X-ray surveys. Based on this, type-I AGN are expected 
to contribute ~ 20 — 30% of the bright 15 fim counts 
(Si5 > 10 mJy, see Fig. ||), and negligible fractions at 
fainter fluxes or longer wavelengths. 

5.5. Contributions of IR galaxies to the ORB 

The 15/im counts in Fig. display a remarkable conver- 
gence below Si5 ^ 0.2 mJy, proven by at least three in- 
dependent surveys. The observed asymptotic slope flatter 
than —2 in differential count units implies a modest contri- 
bution to the integrated CIRB flux by sources fainter than 
0.1 mJy, unless a sharp upturn of the counts would hap- 
pen at much fainter fluxes with a very steep number count 
distribution, which seems rather unplausible. A meaning- 
ful estimate of the CIRB flux can then be obtained from 
direct integration of the observed mid-IR counts: this com- 
putation has been done by Elbaz et al. (2001), who flnd a 
value at 15 //m of 2. 6 ±0.5 nW/m^/sr contributed by LW3 
sources brighter than S15 = 40 //Jy (corresponding to the 
datapoint at 15 /im in Fig. |l|, and consistent with the re- 
sults by Biviano et al. 2001; the other closeby point in the 
figure comes from a similar integration of ISO counts at 
7 /Ltm). From our reference evolutionary model we expect 
that the contribution of fainter sources would bring the 
total background to 3.3 nJy/m^/sr. 

Comparing these values with the upper limits set by 
the observed TeV cosmic opacity (dotted histogram in Fig. 
|l|) confirms that the ISOCAM surveys have resolved a 
significant fraction (50-70%) of the CIRB intensity in the 
mid-IR. 

Unfortunately, we do not have a way to test directly 
how much of the bolometric CIRB these faint 15 /im sur- 
veys contribute. In particular, the depths of the ISO far-IR 
surveys (FIRBACK and FLATS, see Dole et al. 2001) are 
not enough to resolve more than few percent of the CIRB 
at its peak wavelength. 

Using the locally established good correlation between 
the mid-IR and the far-IR fiuxes for local IR galaxies and 
after a careful analysis of the incidence of AGNs among the 



faint ISOCAM sources (this is essential because starbursts 
and AGNs have different IR SEDs, with peak emissions at 
^ 60— 100 /im and '^ 20 — 30 /im, respectively), Elbaz et 
al. (2001) estimate that the sources resolved by CAM LW3 
contribute 60% at least of the COBE/DIRBE background 
at 140 and 240 ^m. 

The good match to the IR multi-wavelength counts 
and related statistics that we found in our previous analy- 
sis by assuming a typical starburst spectrum for the evolv- 
ing population indicates that these IR sources are likely 
dominated by star-formation processes (see also Sect. 6.4). 
Our best-fit model of the multi-wavelength statistics im- 
plies that the ISOCAM sources with S15 > 40 fxiy con- 
tribute a CIRB intensity at 170 /xm of vl{v) ~ 10~^ 
Watt/mVsr, or - 50% of the observed flux (Fig. |). At 90 
/im the fraction rises to ~ 59%. All this supports the con- 
clusion that the population detected by ISO in the mid-IR 
not only contributes a major fraction of CIRB at IBfi, but 
is also likely responsible for a majority contribution to the 
photon energy density contained in the CIRB. 

On the contrary, sources at much higher redshifts, as 
the SCUBA ones are observed to be (see Sect. 6.3 below), 
likely produce a much less significant contribution to the 
bolometric CIRB whatever their comoving volume energy 
production rate might be. The high-z energy production 
can only last for a short cosmic time interval {At oc (1 -|- 
z)~^/^Az), and the generated photons are degraded in 
energy by (1 + z), for a total penalty factor of (1 -I- z)^-^ 
(Harwit 1999). Nevertheless, because of the K-correction, 
high-redshift objects can dominate the CIRB at the longer 
wavelengths, as SCUBA sources are observed to do. All in 
all, the background radiation is not a sensitive tracer of 
the very ancient cosmic phases. 

6. NATURE OF THE FAST EVOLVING IR 
SOURCE POPULATION 

6.1. The mid-IR selected sources 

Given the variety of multi-wavelength imaging data, the 
ISO surveys in the Hubble Deep Fields and the CFRS 
areas (Sect. 3.1) provide ideal tools for tests of the evolving 
population responsible for the upturn of the ISO mid-IR 
counts and for a substantial fraction of the CIRB. 

Aussel et al. (1999 and 2001) report reliably tested 
complete samples of 49 and 63 sources to S15 > 100/zJy in 
the HDF North and South respectively. Flores et al. (1999) 
analyse a sample of 41 sources brighter than S15 ~ 300/iJy 
in the CFRS 1415-1-52 area. These surveys have been per- 
formed with highly redundant spatial and temporal sam- 
pling with ISOCAM LW3, allowing to achieve a precise 
astrometric registering. These sources are identified with 
optical galaxies having typical magnitudes from V 21 
to U ~ 24.5 (Aussel et al. 2001; see also Sect. 3.2). 

The optical-IR SED of a typical faint LW3 source at 
z = 1.14 is reported in Figure |l^, where the LW3 fiux and 
the LW2 upper limit are plotted together with optical- NIR 
fluxes. 
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Fig. 17. Broad-band spectrum of a mid-IR source selected 
by ISOCAM LW3 in the Hubble Deep Field North (Aussel 
et al. 1999), compared with the SED's of M82 (thick con- 
tinuous line), Arp 220 (dashed line), and M51 (dotted 
line). Estimates of the SF rate [based on the M82 and 
Arp 220 templates] and of the stellar mass [based on the 
M51 template] are indicated. 

The dotted line fitting the optical-NIR spectrum and 
corresponding to the SED of a quiescent spiral (M51) falls 
short by a factor ~ 10 of explaining the mid-IR emission, 
whereas SEDs of IR starbursts (Arp220, dashed line; M82 
continuous line) provide more consistent fits to the ob- 
served mid-IR flux after normalizing to the optical/NIR 
spectral intensity. The vast majority of faint ISO sources 
show similar mid-IR flux excesses. 

A clue to the nature of ISO sources can be obtained 
from HST imaging data, providing detailed morphologi- 
cal information, and spectroscopic follow-up available in 
these fields. Flores et al. (1999) find that at least 30 to 
50% of them show evidence of peculiarities and multiple 
structures, in keeping with the local evidence that galaxy 
interactions are the primary trigger of luminous IR star- 
bursts. The Caltech redshift survey in the HDF North by 
Cohen et al. (2000) showed that over 90% of the faint 
LW3 ISO sources are members of galaxy concentrations 
and groups, which they identify as peaks in their redshift 
distributions. It is in such dense galaxy environments with 
low velocity dispersion that interactions produce resonant 
perturbation effects on galaxy dynamics and the most ef- 
ficient trigger of SF. However, both morphological and 
clustering properties of ISO sources need further investi- 
gation, which is currently in progress. 

Flores et al. (1999) report a preliminary analysis of 
optical spectra for IR sources in CFRS 1415-H52, noting 



that a majority of these display both weak emission (OH 
3727) and absorption {Hs) lines, as typical of the e(a) 
galaxy spectral class. Rigopoulou et al. (2000) have ob- 
served with the ISAAC spectrograph on VLT a sample 
of 13 high-z (0.2< z <1.4) galaxies selected in the HDF 
South with Si5 > 100 /iJy: a prominent (EW> 50 A) 
Ha line is detected in almost all of the sources, indicating 
substantial rates of SF after de-reddening corrections, and 
demonstrating that these optically faint but IR luminous 
sources are indeed powered by an ongoing massive dusty 
starburst. 

The e(a) spectral appearence found by Flores et al. 
(1999) is interpreted by Poggianti & Wu (1999) and 
Poggianti, Bressan, Franceschini (2001) as due to selective 
dust attenuation, extinguishing more the newly-formed 
stars than the older ones which have already disrupted 
their parent molecular cloud. These papers independently 
found that 70 — 80% of the energy emitted by young 
stars leaves no traces in the optical spectrum, hence can 
only be accounted for with long- wavelength observations. 

Further efforts of optical-NIR spectroscopic follow-up 
of faint ISO sources are presently ongoing, including at- 
tempts to address the source kinematics and dynamics 
based on line studies with IR spectrographs on large tele- 
scopes (Rigopoulou et al. 2001, in preparation). At the 
moment, for an evaluation of the main properties of the 
IR population we have to rely on indirect estimates ex- 
ploiting the near-IR and far-IR fluxes. We have estimated 
the baryonic mass in stars from fits of template SEDs of 
local galaxies to the observed near-IR broad-band spec- 
trum. Our adopted templates come from the modellistic 
analysis of Silva et al. (1998) of a sample of both inactive 
spirals and starbursts of various masses and luminosities. 
The Initial Mass Function is assumed to be a Salpeter with 
standard low- and high-mass cutoffs (0.15 to 100 Mq). 
Our estimated values of the baryonic masses 10^^ Mq 
at z > 0.4, with 1 dex typical spread, see Fig. |lj) indicate 
that already evolved and massive galaxies preferentially 
host the powerful starbursts. 

For estimating the other fundamental indicator of the 
physical and evolutionary status of the sources - the ongo- 
ing rate of star-formation (SFR) - we have exploited the 
mid-IR flux as an alternative to the (heavily extinguished) 
optical-UV emissions. The capability of the mid-IR flux 
(from both LW3 and LW2 ISOCAM observations) as a 
tracer of the SFR is discussed by Vigroux et al. (1999), 
Elbaz et al. (2001) and Aussel et al. (2001): the fluxes in 
these IR bands appear tightly correlated with the bolo- 
metric (mostly far-IR) emission, which is the most robust 
measure of the number of massive reddened newly-formed 
stars, and is also correlated with the radio emission of stel- 
lar origin. Only very extinguished peculiar sources (e.g. 
Arp 220), for which the mid-IR spectrum is self-absorbed, 
escape this correlation. Note that the ISO mid-IR flux at 
these very faint limits provides advantages over the radio 
to be a more reliable (having a tighter correlation with the 
bolometric flux, see e.g. Cohen et al. 2000) and more sen- 
sitive indicator of star formation. Of the 49 ISO sources 
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Fig. 18. Evaluations of the star formation rates [from an 
estimator based on the mid-IR flux] and baryonic masses 
[from fits of the NIR SED] as a function of redshift for 
galaxies selected by ISOCAM LW3 at IS^im in the HDFN 
(filled squares) and CFRS 1415+52 (open circles). 



in the HDFN only 7 are detected in an ultra-deep radio 
map at 1.4 GHz by Richards et al. (1998): the ISO 15 ^lm 
completeness limit of 0.1 mJy would correspond to a flux 
of few /iJy at 1.4 GHz (i.e. quite below the Richards et 
al. limit of 40 /iJy), taking into account the radio/far-IR 
luminosity correlation for typical starbursts. 

The rates of SF indicated by the fits to the mid-IR 
flux for sources at z > 0.5 range from several tens to few 
hundreds solar masses/yr, i.e. a substantial factor larger 
than found for optically-selected galaxies at similar red- 
shifts (e.g. Ellis 1997). 

Altogether, the galaxy population dominating the faint 
mid-IR counts and substantially contributing to the bolo- 
metric CIRB intensity (assumed typical SB SEDs) ap- 
pears to be composed of luminous (Lboi ~ 10^^ — 10^^-^ 
Lq) starbursts in massive {M ^ 10^^ Mq) galaxies at 
z 0.5 — 1.3, observed during a phase of intense stellar 
formation {SFR ~ 100 Mq/yt). The typically red colors 
of these systems suggest that they are mostly unrelated 
to the faint blue galaxy population dominating the optical 
counts (Ellis 1997), and should be considered as an inde- 



pendent manifestation of (optically hidden) star formation 
(Elbaz 1999, preprint; Aussel 1998). 

6.2. Far-IR selected galaxy samples 

Surveys at longer wavelengths suffer quite severe problems 
in the identification of the optical counterparts of the IR 
sources. 

The FIRBACK/ELAIS surveys have resolved a mod- 
est fraction (- 5%, see Dole et al. 2001) of the CIRB at 
its peak wavelength of 170 /im, the limit being imposed by 
source confusion. Because of the missing information on 
the far-IR LLF, the interpretation of the counts in Figs. 
^, ^ and ^ is subject to some uncertainties. Our best-fit 
multi-wavelength model implies that the observed counts 
at the faint flux limit are a moderate factor (^ 2 — 3) above 
the source areal density corresponding to no-evolution (see 
Fig. ^). The model predicts that the detected sources lye at 
moderate-redshifts (the majority at z < 0.5, see Fig. p2| ). 
The multi-wavelength follow-up performed at 1.4 GHz, 1.3 
mm, 850 and 450 ^m, as well as optical/NIR identifica- 
tions and spectroscopy based on cross-correlations with 
deep radio surveys (e.g. Sanders 2001), seem to show that 
the majority of the sources are local {z < 0.5), with 10% 
or so being found at z ^ 1 or higher (Dole et al. 2001). 

Scott et al. (2000), in particular, have obtained data at 
450 and 850 fj,m for 10 FIRBACK sources with accurate 
radio astrometry: although the combined FIR/radio se- 
lection and the sub-mm follow-up may somehow bias the 
result, the FIR- mm SEDs compared with plausible far-IR 
template spectra tentatively indicate mostly low redshifts 
for these sources, with a minority being at z ~ 1. 

In the future, deeper far-IR observations will be pos- 
sible with SIRTF, while a proper characterization of the 
faint far-IR population will require the Herschel's better 
spatial resolution. 

6.3. Faint millimetric sources 

Thanks to the unique iiT— correction for dusty spectra, 
deep millimetric surveys are capable to detect starburst 
galaxies over an extremely wide redshift interval. However, 
the sensitivities achievable by present-day instruments 
and the modest surveyed areas imply that only sources 
at z > 1 are selected; consequently, only the very high 
luminosity tail of the population of IR sources (essentially 
ULIRGs with Lboi > 10^^ L0) is detectable in this way. 

Therefore, ISOCAM and millimetric telescopes pro- 
vide extremely complementary sampling capabilities in 
terms of redshift coverage (typically z < 1 for sources 
selected by ISO and z > 1 by SCUBA) and source lu- 
minosities (mostly Lboi < 10^^ L0 by ISO, and larger by 
SCUBA). 

Unfortunately, the extreme properties of the mm- 
selected sources entail a dramatic difficulty to identify the 
sources. Several factors contribute: the very high redshifts 
imply that the optical counterparts are extremely faint 
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and red, hence largely unaccessible by optical spectroscopy 
(Small et al. 1999 identify a fraction of SCUBA sources 
with Extremely Red Objects, ERO's, see also discussion 
in Dey et al. 1999). Furthermore, the dominance of dust 
emission implies very extinguished optical-NIR spectra. 
Although the diffraction-limited beams of both SCUBA 
(15 arcsec) and IRAM (11 arcsec) are much sharper than 
the ISO far-IR beam (~60-90 arcsec), the faintness of 
optical counterparts implies a very high chance of mis- 
identification. It is remarkable that, in spite of the impor- 
tant effort for SCUBA source identification, of the 100- 
200 sources only three have reliable identifications and 
redshifts at present. All these three have been observed 
to contain massive reservoirs of gas (Frayer et al. 1999; 
Ivison et al. 2000a). 

Given the extreme difficulty to get the redshift from 
optical spectroscopy, some millimetric estimators have 
been devised. The most promising of such techniques, ex- 
ploiting the S850/j/5'20cm flux ratio as a monotonic function 
of redshift (Carilli & Yun 1999), confirms in a statistical 
sense that faint SCUBA sources are ultra-luminous galax- 
ies at typical z~lto~3 — 4 (Barger et al. 1999; Small 
et al. 2000). The predicted z-distribution by our reference 
model (Fig. ^) is in rough agreement with these esti- 
mates. 

As suggested by several authors (Franceschini et al. 
1994; Lilly et al. 1999; Granato et al. 2001), the similar- 
ity in properties (bolometric luminosities, SEDs) between 
this high-z population and local ultra-luminous IR galax- 
ies argues in favour of the idea that these represent the 
long-sought "primeval galaxies", those in particular origi- 
nating the local massive elliptical and SO galaxies. This is 
also supported by estimates of the volume density of these 
objects in the field ~ 2 — 4 x 10"'' Mpc"'^, high enough 
to allow most of the E/SO to be formed in this way (Lilly 
et al. 1999). As for the E/SO galaxies in clusters, the re- 
cent discovery by SCUBA of a significant excess of very 
luminous {L ~ 10^^ Lq) sources at 850 /.tm around the 
z=3.8 radiogalaxy 4C41.17 (Ivison et al. 2000b) may in- 
dicate the presence of a forming cluster surrounding the 
radiogalaxy, where the SCUBA sources would represent 
the very luminous ongoing starbursts. 

By continuity, the less extreme of the starbursts (those 
with L ~ 10" - 10^2 Lq) discovered by ISOCAM at 
lower redshifts may be related to the origin of lower mass 
spheroids and spheroidal components in later morpholog- 
ical type galaxies. 

6.4. AGN contributions to the source energetics 

A natural question arises as of how much of the bolomet- 
ric flux in these IR/mm sources is contributed by grav- 
itational accretion rather than stellar emission. Almaini, 
Lawrence and Boyle (1999) have suggested that a mini- 
mum of 10 - 20% of the CIRB at 850 /^m (and a sim- 
ilar fraction of the bolometric one) may be due to ob- 
scured AGNs, and that this fraction could even be quite 



larger. Unfortunately, probing the nature of the faint IR- 
mm sources at high-redshifts turned out to be exceed- 
ingly difficult, since the optical-UV-soft-X-ray primary 
emission is almost completely re-processed by dust into 
an IR spectrum barely sensitive to the properties of the 
primary incident one. Ivison et al. (2000a) ffird indications 
for the presence of (type-II) AGN components in three of 
the seven SCUBA sources in cluster fields, a fraction not 
inconsistent with that observed in local ULIRGs of similar 
extreme luminosity. 

Preliminary inspection of the Ha line profiles for the 
faint ISO mid-IR sources (Rigopoulou et al. 2000), to- 
gether with constraints set by the 15 to 7 /zm flux ratios 
(for sources at z ~ 0.5 — 1 this ratio measures the pres- 
ence of very hot dust heated by the AGN, which is absent 
in starbursts), indicate that the large majority of sources 
are also mostly powered by a starburst rather than an 
AGN. Tran et al. (2001) consistently find that the con- 
tribution of gravitational accretion to the IR emission by 
local ULIRGs becomes important, or even dominant, only 
in the very high luminosity regime (Lboi > 10^^'^ L©)- 

An important diagnostic is being offered by observa- 
tions of the hard X-ray flux, since starbursts are weaker 
X-ray emitters than any kinds of AGNs. The Chandra 
X-ray observatory has performed several deep investiga- 
tions of the high-z SCUBA population (Fabian et al. 2000; 
Hornschemeier et al. 2000; and Barger et al. 2001). Only 
a very small percentage of the objects turn out to be 
in common, the two classes of sources being largely or- 
thogonal. The two high-z sub-mm sources detected with 
Chandra by Bautz et al. (2000) were both previously clas- 
sified as AGN based on optical spectra. Unless the large 
majority of sub-mm sources are Conipton-thick and any 
hard X-ray scattered photons are also photo-electrically 
absorbed, the conclusion is that the bulk of the emission 
by high-luminosity SCUBA sources is due to star forma- 
tion (in agreement with a dominant stellar emission in 
local ULIRGs inferred by Genzel et al. 1998). The frac- 
tion of the CIRB at 850 /zm due to AGNs was estimated 
by Barger et al. (2001) to be not likely larger than 10%. 

Fadda et al. (2001) and Elbaz et al. (2001) discuss 
cross-correlations of faint ISO samples with deep hard X- 
ray maps from XMM and Chandra. They find that typi- 
cally 10% of the faint ISO sources show hard X-ray evi- 
dence for the presence of an AGN. 

Although the detailed interplay between starbursts 
and AGNs is still an open issue, it is quite likely that minor 
AGN contributions are present in a substantial fraction of 
the active IR population: Risaliti et al. (2000) and Bassani 
et al (2001 ) find trace AGN emission in 60 - 70% of local 
ULIRGs, based on BeppoSAX hard X-ray data. Similarly, 
a significant fraction of the high-z IR starbursts discovered 
by ISO and SCUBA may contain (energetically-negligible) 
low-luminosity AGNs, detectable in hard X-rays, and pos- 
sibly responsible for the bulk of the X-ray background. 
Still their IR emission is likely to be mostly of stellar ori- 
gin. 
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6.5. Discussion 

Among various samples of faint sources selected at long 
wavelengths, those detected by ISOCAM in deep and 
ultra-deep surveys allow the most precise quantification 
of the cosmic history of the IR population. The ISOCAM 
LW3 extragalactic counts, extending over 4 orders of mag- 
nitude in flux (Fig. ^ when combined with the IRAS 12 
/xm local surveys, provide rather detailed constraints on 
the evolution pattern. The LW3 sources not only con- 
tribute a dominant fraction of the CIRB in the mid- 
IR, but they are also likely important contributors to 
the CIRB at longer wavelengths. The extremely high lu- 
minosities and redshifts and modest volume densities of 
SCUBA sources indicate that they probably produce only 
a small fraction of the bolometric CIRB energy. It should 
not be forgotten, however, that these inferences will re- 
main model dependent untill we will be able to resolve 
into sources a significant fraction of the CIRB at its peak 
wavelengths, and this will have to wait for the operation 
of Herschel/FIRST. 

Because of their non-extreme properties, ISOCAM 
LW3 sources can be fairly unambiguously identified and 
investigated in the optical. The outcome of our spectro- 
scopic observations is that the faint population making up 
the CIRB in the mid-IR appears dominated by actively 
star-forming galaxies with substantial Ha emission. 

The LW3 ISOCAM counts and redshift distributions 
require extremely high rates of evolution of the 15//m lu- 
minosity function up to ~ 1, with preference for evolu- 
tion (in both source luminosity and spatial density) of a 
population of IR starbursts contributing little to the lo- 
cal LF. A natural way to account for this would be to 
assume that this population consists of otherwise normal 
galaxies, but observed during a dust-extinguished lumi- 
nous starburst phase, and that its extreme evolution is 
due to an increased probability with z to observe a galaxy 
during such a transient luminous starburst event. 

The widespread evidence that starbursts are triggered 
by interactions and merging suggests that the number den- 
sity evolution could be interpreted as an increased proba- 
bility of interaction back in time. Assuming that the phe- 
nomenon is dominated by interactions in the field and a 
velocity field constant with z, than this probability would 
scale roughly as oc p{z)'^ cx (1 -f z)^, p being the number 
density in the proper (physical) volume. A more complex 
situation is likely to occur, since also the velocity field 
evolves with z in realistic cosmological scenarios and if we 
consider that the most favourable environment for inter- 
actions are galaxy groups, which are observed to host the 
majority of ISOCAM distant sources (Cohen et al. 2000). 
In turn, the increased luminosity with z of the typical 
starburst is due, qualitatively, to the larger amount of gas 
available in the past to make stars. 

How this picture of a 2-phase evolution of faint IR 
sources compares with results of optical and ncar-IR deep 
galaxy surveys has not been investigated by us. Since, be- 
cause of dust, most of the bolometric emission during a 



starburst comes out in the far-IR, we would not expect 
the optical surveys to see much of this violent starbursting 
phase revealed by IR observations. Indeed, B-band counts 
of galaxies and spectroscopic surveys are interpreted in 
terms of number-density evolution, consequence of merg- 
ing, and essentially no evolution in luminosity. The Faint 
Blue Object population found in optical surveys can be in- 
terpreted in our scheme as the post-starburst population, 
objects either observed after the major event of SF, or 
more likely ones in which the moderately extinguished in- 
termediate age (> 10^ yrs) stars in a prolonged starburst 
(several 10^ yrs) dominate the optical spectrum. In this 
sense optical and far-IR selections trace different phases 
of the evolution of galaxies, and provide independent sam- 
pling of the cosmic star formation. 

A lively debate is currently taking place about the ca- 
pabilities of UV-optical observations to map by themselves 
the past and present star-formation, based on suitable cor- 
rections for dust extinction in distant galaxies. Adelberger 
et al. (2000) suggest that the observed 850 pm galaxy 
counts and the background could be explained with the 
optical Lyman drop-out high-z population by applying a 
proportionality correction to the optical flux and by taking 
into account the locally observed distribution of mm-to- 
optical flux ratios. 

On the other hand, a variety of facts indicate 
that optically-selected and IR/mm-selected faint high- 
redshift sources form almost completely disjoint samples. 
Chapman et al. (2000) observed with SCUBA a subset of 
z ~ 3 Lyman-break galaxies having the highest predicted 
rates of SF as inferred from the optical spectrum, but de- 
tected only one object out of ten. van der Werf et al. (2000) 
found that the procedures adopted to correct the optical- 
UV spectrum and to infer their sub-mm fluxes failed in 
the case of two strongly lensed Lyman-break galaxies ob- 
served with SCUBA. A similar dichotomy is observed in 
the local universe, where the bolometric flux by luminous 
IR galaxies is mostly unrelated with the optical emission 
spectrum (Sanders & Mirabel 1996). 

Rigopoulou et al. (2000), Poggianti & Wu (2000) and 
Poggianti et al. (2001) report independent evidence on 
both local and high-z luminous starbursts that typically 
70% to 80% of the bolometric Hux from young stars leaves 
no traces in the UV-optical spectrum, because it is com- 
pletely obscured by dust. As there seems to be no robust 
" a priory" way to correct for this missing energy, we con- 
clude that only long-wavelength observations, with the ap- 
propriate instrumentation, can eventually measure SF in 
galaxies at any redshifts. 



7. GLOBAL PROPERTIES: THE COMOVING 
SFR DENSITY AND CONSTRAINTS FROM 
BACKGROUND OBSERVATIONS 
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Fig. 19. Evolution of the comoving luminosity density 
for the IR-selected population based on the model of IR 
evolution discussed in Sect 5. The luminosity density is 
expressed here in terms of the star formation rate den- 
sity (computed from the far-IR luminosity assuming a 
Salpeter IMF, according to the recipes reported in Rowan- 
Robinson et al. (1997, their eq. [7]). The IR evolution is 
compared with data coming from optical observations by 
Lilly et al. (1996), Connolly et al. (1997) and Madau et al. 
(1996), transformed to our adopted Qy^ = 0.3, = 0.7 
cosmology. Dotted line: quiescent population. Short-dash 
line: evolving starbursts. Long dashes: type-I AGNs (with 
arbitrary normalization). The continuous line is the to- 
tal. The shaded horizontal region is an evaluation of 
the average SFR in spheroidal galaxies by Mushotzky & 
Loewenstein (1997). 



30% of the bolometric emission by young stars in IR galax- 
ies, see Sect. 6). 

The fast evolution inferred from the IR observation 
should however level off at z > 1 , to allow consistency with 
the observed z-distributions for faint ISOCAM sources 
(Fig. se also Chary & Elbaz 2001). Another important 
constraint in this sense comes from the observed spectral 
shape of the CIRB, with its apparent sharp peak around 
A ~ 100 — 200 iim and the fast turnover longwards: fitting 
it with standard dust-emission spectra implies a maximum 
in the comoving galaxy IR emissivity close to z = 1. This 
is also confirmed by attempts to derive the average time- 
dependent IR volume emissivity from deconvolution of the 
CIRB spectrum, see Gisper et al. (2000) and Takeuchi et 
al. (2001). 

Surveys of Hq, line emission (the best optical indicator 
of SF) from high-z galaxies indicate similar evolution of 
the SFR density, with a similarly sharp change in slope 
occurring at z ~ 1 (see van der Werf, Moorwood and Yan 
2001 for a review). Our present results offer the advantage, 
however, to be unaffected by the uncertain extinction cor- 
rections of the optical SF indicators, and to exploit the 
very robust constraint set by the CIRB. 

Altogether these results indicate that the history of 
galaxy long-wavelength emission does probably follow a 
general path not much dissimilar from that revealed by 
optical-UV observations, by showing a similar peak activ- 
ity around 2 ~ 1 — 1.5, rather than being confined to the 
very high-redshifts as sometimes was suggested based on 
SCUBA results. This confirms that the bulk of the galaxy 
activity is to be placed around z = 1, which is evident from 
Fig. |l^ if the dependence of the cosmological timescale on 
redshift is considered (see also Harwit 1999 and Haarsma 
& Partridge 1998). 

As a final note, the estimated rate of evolution of the 
IR volume emissivity of galaxies appears in Fig. |l^ to be 
even higher than the evolution rate for type-I AGNs. 



7.1. Evolution of the comoving luminosity density and 
SFR 

We are not presently in the position to derive an indepen- 
dent assessment of the evolutionary SFR density based on 
the available complete samples of faint IR sources, since 
the process of optical identifications is far from complete. 

Again, only model-dependent estimates of the SFR 
density as a function of redshift are possible at the mo- 
ment. The prediction based on our previously described 
evolution scheme is reported in Fig. There is a clear in- 
dication here that the contribution of IR-selected sources 
to the luminosity density significantly exceeds those based 
on optically selected sources, and that the excess may be 
progressive with redshift up to z ~ 1. Consider tliat, in 
any case, the optical and IR estimates of the SFR in Fig. 
|T| refer to largely distinct source populations (extinction- 
corrected optical fluxes account on average for only 20- 



7.2. Energy constraints from background observations 

Further constraints on the high-redshift far-IR/sub-mm 
population can be inferred from observations of the global 
energetics residing in the CIRB and optical backgrounds. 
These imply a very substantal demand on contributing 
sources, as detailed below in schematic terms. 

Let us assume that a fraction of the universal 

mass density in baryons pb — ^—^s^^-^^h — 7 10^*^(1 -I- 
z)^hlQQh [Mo/Mpc'^] undergoes at redshift z* a transfor- 
mation (either processed in stars or by gravitational fields) 
with radiative efficiency e. The locally observed energy 
density of the remnant photons is 



P7 = Pb 
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rib /* 2.5 
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For stellar processes, e is essentially determined by the 
IMF: e = 0.001 for a Salpeter IMF and a low-mass cutoff 
M„ii„ = 0.1 Mq, e = 0.002 and e = 0.003 for Af^in = 2 
and Mmin = 3, while e gets the usually quoted value of 
e = 0.007 only for M^in = 10 Mq. 

7.2.1. Constraints from the integrated optical 
background 

Let us adopt for the optical/near-IR bolometric emission 
by distant galaxies between 0.1 and 7 /xm the value given 
in eq.(^). We discussed evidence that in luminous star- 
bursts the optical spectra are only moderately contributed 
by the starburst emission itself, the latter being largely 
hidden in the far-IR. Then let us assume that the opti- 
cal/NIR background mostly originates by moderately ac- 
tive SF in spiral disks and by intermediate and low-mass 
stars. As observed in the Solar Neighborhood, a good ap- 
proximation to the IMF in such relatively quiescent en- 
vironments is the Salpeter law with standard low-mass 
cutoff, corresponding to a mass-energy conversion efff- 
ciency e 0.001. With these parameter values, we can 
reproduce the whole optical BKG intensity of eq.(^ by 
transforming a fraction ~ 10% of all nucleosynthetic 
baryons into (mostly low-mass) stars, assumed the bulk of 
this process happened at ~ 1.5 and 5% of the closure 
value in baryons (for our adopted Hq = 50 Km/ s/Mpc, or 
r^b/iioo — 0.012, consistent with the theory of primordial 
nucleosynthesis) : 

i^^('^)|opt ^ 

20 10-9/i2 ilL A ( ^Watt/mVsr. 
^"0.05 0.1 Vl + 2, y 0.001 ' ' 

A local density in low-mass stars is generated in this way 
consistent with the observations (Ellis et al. 1996), and 
based on standard mass to light ratios: 

Pb(stars) ~ 7 lO^^f^Qb - 3.4 10^ Mo/Mpc^ (12) 

For typical solar metallicities, this corresponds to a local 
density in metals of 

pz (stars) 

1.67 10V*i|-^^b Mq/Mpc^ ~ 8 lO^Me/Mpc^. (13) 

7.2.2. Explaining the CIRB background 

Following our previous assumption that luminous star- 
bursting galaxies emit negligible energy in the optical- 
UV and most of it in the far-IR, we coherently assume 
that the energy resident in the CIRB background (eq. |^ ) 
originates from star- forming galaxies at median ~ 1.5. 
The amount of baryons processed in this phase and the 
conversion efficiency e have to account for the combined 
constraints set by eqs.(|l^), eqs.(|3|) and (||), that is to pro- 
vide a huge amount of energy without contributing much 



stellar remnants to the locally observed amount. A viable 
solution is then to change the assumptions about the stel- 
lar IMF characterizing the starburst phase, for example to 
a Salpeter distribution cutoff below M^in = 2 Mq, with a 
correspondingly higher efficiency e = 0.002 (see Sect. 7.2). 
This may explain the energy density in the CIRB (eq. |^ ) : 

i^^(j^)|fir ^ 

40 10-9/i2 Al A ( AA^ Watt/mVsr, 
^° 0.05 0.1 \l + zj 0.002 ' ' ' 

assumed that a similar amount of baryons, /* ~ 10%, as 
processed with low efficiency during the "inactive" secular 
evolution, are processed with higher radiative efficiency 
during the starbursting phases, producing a two times 
larger amount of metals: p{metals) ^ 1.5 10^ Mq/Mpc'. 
Note that by decreasing M,„in during the SB phase would 
decrease the efficiency e and increase the amount of pro- 
cessed baryons /*, and would bring to exceed the locally 
observed mass in stellar remnants (eq.[Q). 

The above scheme is made intentionally extreme, to 
illustrate the point. The reality is obviously more com- 
plex, e.g. by including a flattening at low mass values in 
the Salpeter law (e.g. Zoccali et al 1999) for the solar- 
neighborhood SF and, likewise, a more gentle convergence 
of the starburst IMF than a simple low-mass cutoff. 

7.3. Galactic winds and metal pollution of the 
inter-cluster medium 

A direct prediction of the above scheme is that most of 
the metals produced during the starburst phase have to 
be removed by the galaxies to avoid exceeding the locally 
observed metals in galactic stars (eq.||l^). There is clear 
evidence in local starbursts, based on optical and X-ray 
observations, for large-scale super-winds out-gassing high- 
temperature enriched plasmas from the galaxy. Our expec- 
tation would be that the amounts of metals originating 
from the SF processes producing the CIRB are hidden in 
hot cosmic media. 

Where all these metals are? Most likely the polluted 
plasmas are hidden in the diffuse (mostly primordial and 
un-processed) inter-cluster medium with densities and 
temperatures preventing to detect them. On the other 
hand, an interesting support to our scheme is provided 
by observations of rich clusters of galaxies, considered as 
closed boxes from a chemical point of view, as well as rep- 
resentative samples of the universe. The mass of metals in 
the ICP plasma can be evaluated from the total amount 
of ICP baryons (~5 times the mass in galactic stars) and 
from their average metallicity, ~40% solar. The mass of 
ICP metals is then Mmetaisjcp ^ 5 x 0.4 {Z/Zq) M^tars, 
i.e. two times larger than the mass of the metals present 
in galactic stars and consistent with the mass in metals 
we inferred to be produced during the SB phase. 

Hence, the same starbursts producing the ICP metals 
are also likely responsible for the origin of the CIRB. In 
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a similar fashion, Mushotzky & Loewenstein (1997) used 
their metaUicity measurements of clusters to estimate the 
contribution of spheroidal galaxies to the SFR density (see 
Fig.[l|). 

7.4. A two-phase star-formation: origin of galactic 
disks and spheroids 

Our reference model implies that star formation in galax- 
ies has proceeded along two phases: a quiescent one tak- 
ing place during most of the Hubble time, slowly building 
stars with standard IMF from the regular flow of gas in 
rotational supported disks; and a transient actively star- 
bursting phase, recurrently triggered by galaxy mergers 
and interactions. During the merger, violent relaxation 
redistributes old stars, producing de Vaucouleur profiles 
typical of galaxy spheroids, while young stars are gener- 
ated following a top-heavy IMF. 

Because of the geometric (thin disk) configuration of 
the diffuse ISM and the modest incidence of dusty molec- 
ular clouds, the quiescent phase is only moderately af- 
fected by dust extinction, and naturally produces most of 
the optical/NIR background (including NIR emission by 
early- type galaxies completely deprived of an ISM). 

The merger-triggered active starburst phase is instead 
characterized by a large-scale redistribution of the dusty 
ISM, with bar-modes and shocks, compressing a large frac- 
tion of the gas into the inner galactic regions and trigger- 
ing formation of molecular clouds. As a consequence, this 
phase is expected to be heavily extinguished and the bulk 
of the emission to happen at long wavelengths, naturally 
originating the cosmic IR background. Based on dynami- 
cal considerations, we expect that during this violent SB 
phase the elliptical and SO galaxies are formed in the most 
luminous IR SBs (corresponding to the SCUBA source 
population), whereas galactic bulges in later- type galax- 
ies likely originate in lower IR luminosity starbursts (the 
ISO mid-IR population). 

The presently available IR-selected galaxy samples, 
dominated as they are by K-correction and selection ef- 
fects, cannot allow to establish the precise evolutionary 
timescales as a function of source luminosity. In our best- 
fit model, both SCUBA-selected ULIRG and ISO-selected 
LIRG galaxies have the same evolution history: if any- 
thing, SCUBA sources originating massive E/SOs might 
evolve on a faster cosmic timescale. This could be still in 
line with the expectations of hierarchical clustering mod- 
els if we consider that SCUBA sources likely trace the very 
high-density environments (galaxy clusters) with an accel- 
erated merging rate at high-z, while ISO sources are likely 
related with lower-density environments (galaxy groups or 
the field) entering the non-linear collapse phase at later 
cosmic epochs (e.g. Franceschini et al. 1999). 

Finally, if indeed the IMF characteristic of the star- 
burst phase is deprived of low-mass stars, as suggested 
in the previous paragraphs, a consequence would be that 
the excess blue stars formed during the starburst would 



quickly die and disappear, leaving the colors of the emerg- 
ing remnant as typically observed for early-type galaxies 
and keeping consistent with the evidence that the stel- 
lar mass in spheroidal galaxies does not change much for 
z<l. 

8. SUMMARY AND CONCLUSIONS 

We have analyzed a large dataset derived from deep galaxy 
surveys at long wavelengths, exploiting in particular new 
ISO and published SCUBA observations, but also includ- 
ing data from IRAS, COBE and IRAM. This study of 
galaxy evolution at long wavelengths benefit also by the 
unique situation to combine constraints coming from faint 
resolved sources with data on the integrated source emis- 
sion provided by the spectral intensity of the cosmic IR 
background. The main results of our analysis are hereby 
summarized. 

- (1) The deep surveys by ISOCAM LW3 at 15 ^m pro- 
vide the most precise quantification of statistical prop- 
erties of faint IR galaxies, in the form of deep source 
counts and redshift distributions for faint complete 
samples. The numerous ISO-selected 15 /im sources 
also allow to perform detailed physical investigations 
of the high-z IR population, since their optical counter- 
parts are rather straightforward to identify. Analyses 
of the mid-IR to far-IR flux correlations suggest that 
not only at the ISOCAM LW3 flux limit a major frac- 
tion of the mid-IR CIRB intensity is resolved into 
sources, but also that the same sources are likely to 
contribute a substantial fraction of the bolometric en- 
ergy density in the CIRB (Elbaz et al. 2001). Our 
present comparative study of the multi-wavelength 
counts and LF's entirely confirms this conclusions, by 
showing that the various statistics on IR sources can 
only be reconciled each other by assuming IR SEDs 
typical of starbursts, in which case ~ 50% of the bolo- 
metric CIRB is contributed by ISOCAM 15 /Ltm sources 
brighter than S15 = 40 /iJy. The sources of this impor- 
tant cosniological component, including ~ 70% of the 
integrated bolometric emission by galaxies, can then 
be investigated in the ISOCAM population. By con- 
trast, our model predictions suggest that the extreme 
luminosities and redshifts and low spatial density of 
the SCUBA-selected population are such that only a 
modest fraction of the bolometric CIRB energy can 
be produce by them, although they can dominate it at 
the long- A. For a constant comoving volume emissivity 
as a function of z (see Fig. ^9|), high-redshift sources 
suffer a (1 -|- z)^^'^ penalty factor in their contribution 
to the CIRB. 

— (2) The most robust conclusion coming from the 
present analysis appears to be that of a very rapid 
increase of galaxy long- wavelength volume emissivity 
with redshift, paralleled by an increased incidence in 
high-redshift sources of dust extinction and thermal 
dust reprocessing, with respect to locally observed 
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galaxies. This is the faster evolution rate observed 
for galaxies at any wavelengths [p{L,z) oc {1 + z)^ 
if averaged over the whole galaxy population, see Fig. 

and is higher than the rates inferred for quasars. 
We have found that the shape of the LW3 counts (a 
roughly Euclidean behavior down to S15 ~ 5 niJy fol- 
lowed by a very fast upturn) is better fit by assuming 
that for only a fraction of local galaxies the IR emissiv- 
ity evolves back in time. In this case strong evolution 
both in number density and in the average source lu- 
minosity is indicated. 

(3) The combined constraints set by the z-distributions 
and by the spectral shape of the CIRB impose that this 
fast evolution rolls-over around rcdshift 1, and keeps 
flat above. The evolution of the luminosity density is 
so rapid up to z ~ 1 that no much further increase is 
left for the higher redshifts. Consequently, the history 
of galaxy formation traced by long-wavelength obser- 
vations does not appear fundamentally different from 
the one inferred from optical observations, and shows 
a maximum around z ^ 1. Scenarios in which a sub- 
stantial fraction of stellar formation happens at very 
high-z (e.g. producing the bulk of stars in spheroidal 
galaxies at z > 2—3) are not supported by our analysis, 
and appear to conflict in particular with the observed 
shape of the CIRB at 1000 < A < 100 ^m. Obviously, 
for a deflnitive proof of these statements it will be nec- 
essary to resolve a signiflcant fraction of the CIRB at 
its peak (A ~ 100 to 300 /im). The next important 
step in this sense is expected from the Herschel Space 
Observatory. 

(4) Our present results, together with preliminary 
spectroscopic studies of the faint IR sources, suggest 
that only a minor fraction of their IR flux originates 
from AGN activity, the bulk of it being likely due to 
star formation. Our best-fits to the SEDs indicate mas- 
sive systems hosting violent starbursts {SFR ^ 100 
M0/yr). We have suggested that the most natural in- 
terpretation for the strong observed evolution is to as- 
sume that the evolving starbursting population con- 
sists of otherwise normal galaxies observed during a 
dust-extinguished short-lived but luminous starburst 
event. The strong increases with redshift of the prob- 
ability of interactions (as partly due to a plain geo- 
metrical effect in the expanding universe) and of the 
effects of interactions (due to the more abundant fuel 
available in the past), likely explain the observed rapid 
evolution. Galaxy interactions and mergers are empha- 
sized by the present analysis as a crucial cosmogonic 
driver for galaxy evolution. 

(5) Our suggested evolutionary scheme considers bi- 
modal star-formation in galaxies: SF during a long- 
lived quiescent phase, and enhanced SF taking place 
in transient starburst phases recurrently triggered by 
interactions and merging. The former would be respon- 
sible for building of galaxy disks, the latter for the 
assembly of spheroidal components in galaxies. Very 
schematically, we have attributed the optical back- 



ground to emission during the quiescent mode, and 
explain the CIRB as mostly due to dust-extinguished 
emission by young stars during the starbursting phase. 

— (6) The large energy content in the CIRB is not easy 
to explain, if we consider the modest efficiency of the 
mass-energy transformations allowed by stellar evolu- 
tion. The most obvious way to alleviate the combined 
constraints set by the local observed amounts of low- 
mass stars and metals would be to assume that a stellar 
IMF somewhat deprived in low-mass stars is character- 
istic of the IR starburst phase. In any case, we expect 
that as much as 2 times the amount of metals in stars 
should be present in the diffuse intergalactic medium 
as a remnant of the ancient SB phase. A support to this 
concept is provided by observations of the abundance 
of metals in hot Intra-Cluster plasmas. 

— (7) Many of the phenomena revealed by long- 
wavelength observations were largely unexpected 
based on UV-optical-NIR observation: among others, 
the rate of cosmic evolution of IR galaxies, the CIRB 
energetics, the luminosities and rates of SF for LIRGs 
and ULIRGs. The attempts so far to infer the IR prop- 
erties based on UV-optical observations are producing 
modest results: there does not seem to be an alter- 
native to long-wavelength observations if we aim at an 
exhaustive and reliable description of the history of SF 
in galaxies. 
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APPENDIX 
PREDICTIONS FOR FUTURE TESTS 

We report in this Section predictions based on the ref- 
erence model, which could be useful for testing it, as well 
as for planning of future observations. 

The confusion noise, which is the fundamental limit- 
ing factor for space IR instrumentation, is evaluated in 
this Section as well as in the paper from the criterion to 
accept sources down to a flux limit corresponding to 1/27 
independent beams. For Euclidean counts, this roughly 
corresponds to a 3f7 confidence limit. Note however that 
for non-Euclidean counts, a case often encountered in the 
infrared, this criterion may under- or over-predict the 3a 
limit according to the count slope (Franceschini 2000). 

SIRTF will operate an imaging camera (MIPS) with 
broad-band filters centered at 24, 70 and 170 /.tm. Fig. |2^ 
(see also Fig. ^ reports predicted counts at these wave- 
lengths. 

The Herschel Observatory (former ESA Cornerstone 
FIRST) will characterize the far-IR (70 to 500 /xm) emis- 
sion by galaxies at any redshifts, by pushing down in flux 
the limit of confusion with its large 3.6m primary mirror. 
Fig. ^ illustrates galaxy counts and the confusion thresh- 
old in two of the Herschel long-A channel at 250 and 450 
fim (other information can be retrived from Figs. |^ and 

The distributions of redshifts for flux-limited sam- 
ples with complete identification are the other funda- 
mental statistical observable for faint distant sources. 
Identifications of FIRBACK/ELAIS sources detected by 
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Fig. 20. Integral number counts at AeS = 24 and 70 
/zm. The SIRTF confusion limits are referred to as the 85 
cm lines. For comparison, an 8m New Generation Space 
Telescope observing at 24 fxin would be confusion limited 
at a source areal density of 10^ sources /sq.deg. (at the 
level of the top axis). 



ISO at 170 /im will require an extensive effort of deep 
radio imaging to reduce the errorbox. A similar effort 
will be required to identify SCUBA- or IRAM-detected 
high-z galaxies. Alternatively, the errorbox will be reduced 
by following-up SCUBA sources with mm interferometers 
(Plateaux de Bure, ALMA). We report in Fig. ^ pre- 
dicted z-distributions for both kind of surveys at the typ- 
ical limiting fluxes. 

Finally, we report in Figures ^ z-distributions for 
confusion-limited surveys by SIRTF and Herschel at 24 
and 450 /im, respectrively. In both cases, the K-correction 
for typical starburst spectra plays in favour of the detec- 
tion of galaxies well above z = 1. Note, in particular, the 
secondary peak at 2: ~ 2 in D{z) for the 24 /im selection. 



Fig. 21. Integral number counts at X^g = 250 and 450 /im, 
relevant for future surveys with the Hershel Observatory. 
See also caption to Fig. 
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due to the PAH emission bundle in the rest-frame ^ 8 /xm 
entering the observable bandwidth at such redshift, pro- 
vides an attractive feature of the forthcoming SIRTF sur- 
veys (e.g. Lonsdale 2001; Lonsdale et al. 2001; Dickinson 
2001). 




Fig. 22. Distributions of redshifts for flux-limited samples Fig. 23. Distributions of redshifts for flux-limited samples 
at 170 (top panel) and 850 fim (bottom), at the limits of at 24 (top panel) and 450 fj,m (bottom), at the limits of 
confusion for ISO and SCUBA. confusion for SIRTF and Herschel. 



